﻿ Gheorghe RUSU Marius SPÎNU Ioana CIOTIR MATEMATICI APLICATE ÎN ECONOMIE IAŞI, 2019 1 Introducere Disciplina „Matematici aplicate în economie” este una dintre disciplinele generale pe care studenții le studiază pentru ca pregătirea lor să fie completă. Ea se alătură altor discipline cu același caracter (limbile străine, dreptul, economia etc.) fără de care un specialist în economie nu poate să își exercite profesia. Trebuie precizat că în alte țări din Europa sau din afara continentului nostru, volumul de cunoștințe matematice predat studenților este mult mai mare și ocupă programe de studii pe mai multe semestre. Ceea ce a fost introdus în această programă și s-a concretizat în prezentul curs reprezintă un minim și un strict necesar. Disciplina este pe de altă parte una obligatorie. Ea este una dintre acele discipline cu credite și fără acestea promovarea nu este posibilă. Numărul de credite în cazul disciplinei de față este egal cu cel al celorlalte discipline obligatorii (adică 5 credite). Este o nouă și importantă dovadă a locului care i se acordă în planul de învățământ al facultății noastre. Ar fi ideal dacă toți studenții noștrii ar avea cunoștințe de matematici elementare echivalente cu cele ale unui absolvent de liceu cu profil economic. Cum acest lucru nu este posibil, mulți studenți provenind din licee cu alte profile unde matematica are o pondere mai mică, ne-am străduit să facem cât mai puține referiri la cunoștințele anterioare de matematică ale studenților. În consecință am cuprins în manual aproape tot ceea ce este necesar pentru parcurgerea și înțelegerea textului. Elementele care sunt preluate pot fi aprofundate și clarificate în activitățile tutoriale. Obiectivele cursului sunt multiple. Pe de o parte ne propunem să familiarizăm studenții cu modalitățile uzuale de modelare matematică a unor fenomene economice. În acest scop vom obișnui studenții să extragă elementele caracteristice esențiale ale unui fenomen economic și să le transpună în relații matematice. Acestei prime etape de modelate îi va urma etapa studiului modelului matematic rezultat și, în final, interpretarea economică a rezultatelor. Pe de altă parte vom face cunoscute studenților două mijloace puternice de optimizare folosite în economie (metode algebrice și metode bazate pe derivată) care au permis rezolvarea unor probleme clasice ale economiei (problema gestionării resurselor, problema transportului, problema repartizarii forței de muncă etc.). 2 Nu în ultimul rând, studenților li se formează obișnuința de a gândi și a acționa respectând regulile stricte ale matematicii. Nu ne propunem să împingem la extrem cerințele privind exactitatea calculului, dar nu trebuie să uităm că avem de-a face cu viitori economiști care vor utiliza mijloace materiale și vor purta răspunderea pentru eficiența și corectitudinea muncii lor. Coroborând cunoștințele de matematică cu cele de informatică, absolvenții noștri vor putea crea aplicații pentru computer care vor completa soft-ul deja existent în domeniul economic. Ei vor deveni nu doar utilizatori ci și creatori de soft economic, ceea ce reprezintă o treaptă calitativă net superioară. Materialul prezentat este structurat în nouă unități de studiu. Această împărțire, împreuna cu planul calendaristic al activităților tutoriale, permite gruparea în trei părți, fiecăreia corespunzându-i o lucrare de evaluare. Conținutul acestor lucrări va privi aspectul teoretic tratat în curs, dar și întrebări de control sub formă de grilă la care studenții vor trebui să răspundă. Întrebările de control la autoevaluare se află în materialul existent pe Blackboard și este accesibil eșalonat studenților, etapă cu etapă. Pentru verificarea cunoștințelor și evaluarea acestora s-au prevăzut trei teme de control. Prima temă de control privește elementele de algebră pe care studenții le parcurg în primele trei unități de studiu. Ea se intitulează astfel „Elemente de algebră liniară”. A doua temă de control privește metodele algebrice de obținere a soluțiilor optime ale problemelor de economie. Această temă se intitulează “Programare liniară și algoritmi asociați” și utilizează material parcursă în unitățile de studio numerele 4, 5 și 6. Ultima temă de control privește material parcursă în unitățile de studiu numerele 7, 8 și 9 și este intitulată ”Elemente de analiză matematică aplicate în economie”. Studenții vor primi temele de control în cadrul activităților tutoriale și le vor depune la termenele indicate la secretarialul ID. Pentru a asigura o asimilare ușoară a materialului teoretic, cursul și activitățile tutoriale cuprind nenumărate exemple. Nu este de dorit și nu încurajăm memorarea fără aprofundare a unor rezultate teoretice. 3 Nu sunt necesare mijloace auxiliare (calculatoare de buzunar, tabele de valori, indrumare de calcul etc.), ci numai materialul editat. Evaluarea se face pe baza temelor de control (60% și a evaluării finale 40%) cu condiția ca ambele să aibă măcar nota 5 (cinci). Dacă notăm cu M media celor trei lucrari de control și cu EvF nota la evaluarea finală (M ≥ 5, EvF ≥ 5), atunci nota finală este . 4 Cuprins 1. Introducere..................................................................................................................... 2 2. Cuprins......................................................................................................................... 5 3. US1 Transformări elementare....................................................................................... 6 4. US2 Spații liniare............................................................................................................ 20 5. US3 Forme pătratice....................................................................................................... 35 6. US4 Programare liniară. Considerații teoretice................................................................... 40 7. US5 Metoda Simplex...................................................................................................... 50 8. US6 Probleme de transport.......................................................................................... 65 9. US7 Serii de numere reale........................................................................................... 80 10. US8 Serii de puteri. Dezvoltarea unei funcții în serie de puteri................................... 107 11. US9 Funcții de mai multe variabile............................................................................. 115 12. Bibliografie................................................................................................................... 134 5 US1 Transformări elementare Nenumărate probleme din matematică utilizează noțiunea de matrice iar transformările elementare de linii aplicate acestora pot înlocui calculele lungi și chinuitoare pe care unele metode algebrice le implică. Ne propunem în această unitate de studiu să prezentăm transformările elementare de linii într-o manieră practică și unitară, și insistăm asupra aplicațiilor lor în determinarea rangului, în calculul determinantului și a inversei unei matrici. De asemenea aplicăm transformările elementare pentru discuția compatibilității și determinarea soluțiilor sistemelor de ecuații liniare. Putem astfel să evităm rezultate clasice (teoremele lui Rouché, Kronecker – Capelli și Cramer) și metode complicate care implică prea multe calcule. Studentul trebuie să se obișnuiască mai întâi să lucreze cu liniile unei matrici, după care trebuie să își stabiliească obiectivul: formarea pe coloana „i” a coloanei „k” din matricea unitate . În continuare trebuie stabilite operațiile necesare scopului propus și efectuarea calculelor implicate de acestea. Ulterior, studentul aplică efectiv metoda însușită pentru determinarea rangului sau inverse unei matrice sau pentru discuția și rezolvarea unui sistem de ecuații algebrice liniare. Timpul necesar însușirii materialului conținut este de aproximativ 6 ore (2 ore pentru învățarea metodei, 2 ore pentru determinarea rangului și a inversei, și două ore pentru aplicații la discuția sistemelor). Metoda este prezentată pe larg în materialul nostru dar și în lista bibliografică la pozițiile [3], [4], [5] și [6]. * * * Studiul matricelor este legat de discuţia şi rezolvarea sistemelor de ecuaţii algebrice liniare şi constituie elementul de control al aproape întregii programe analitice a disciplinei „Matematici aplicate în economie”. Noutatea pe care o propunem se numeşte transformare elementară efectuată asupra liniilor unei matrice, obţinând din aceasta o matrice echivalentă din punctul de vedere al rangului. Departe de a fi o abstractizare inutilă, folosirea transformărilor 6 elementare poate duce, fără calcule dificile, la determinarea rangului, la obţinerea inversei, sau, la obţinerea soluţiei unui sistem algebric liniar. Se cunosc dificultăţile pe care le întâmpinăm la discuţia sau rezolvarea acestui tip de sisteme. Metodele bazate pe teoremele lui Cramer, Kroneker-Capelli sau Rouché sunt greoaie şi folosesc determinanţii, al căror calcul este laborios. Prin transformarea matricei lărgite a unui astfel de sistem se obţine forma explicită din care se pot trage concluziile privind tipul şi soluţiile sale. O categorie specială de soluţii ale unui sistem liniar compatibil nedeterminat o formează soluţiile de bază care vor juca un rol deosebit în rezolvarea problemelor de programare liniară. 7 1.1 MATRICE 1.1.1 RANGUL UNEI MATRICE Prin matrice înţelegem o aplicaţie A:I x J , unde I 1,2,...,m ; J 1,2,...,n , o mulţime oarecare. Ea poate fi reprezentată printr-un tablou de elemente din , aşezate pe linii şi coloane. Pentru o matrice de tipul m x n vom folosi notaţia: 11 12 21 a a 22 a a A ..... ..... am1 am2 ..... ..... ..... ..... 1n 2n a a ..... (1.1.1) amn Suprimând din matricea A, m – r linii şi n – r coloane, obţinem o matrice pătratică de ordinul r. Determinantul acestei matrice se numeşte minor de ordinul r al matricei. Definiţia 1.1.1. Numim rangul unei matrice A de tipul m x n (m n ) un număr natural r cu proprietăţile: a) în matrice există cel puţin un minor de ordinul r diferit de zero. b) toţi minorii de ordinul r + 1 sunt nuli. Observaţie. Se demonstrează că dacă toţi minorii de ordinul r + 1 sunt nuli, atunci sunt nuli şi toţi minorii de ordin mai mare ca r + 1. În continuare vom pune în evidenţă operaţii pe care le vom numi transformări elementare, care se bucură de proprietatea că nu modifică rangul unei matrice. Ele au la bază proprietăţile determinanţilor. 1.1.2 TRANSFORMĂRI ELEMENTARE Definiţia 1.1.2. Numim transformare elementară aplicată unei matrice una din următoarele operaţii: (T1) – înmulţirea unei linii cu un scalar real nenul ( 0); (T2) – adunarea unei linii la o altă linie; 8 (T3) – schimbarea a două linii între ele. Teorema 1.1.1. Transformările elementare nu modifică rangul unei matrice. Demonstraţie. Fie matricea A de tipul m x n şi matricele T(i 1,3) care se obţin din i matricea unitate Im. Înmulţind la stânga matricea A cu matricele T obţinem matricele A T A. i i i Presupunem că rangul matricei A este rA . Vom avea: r i min (r i ,r ) = min (m, r ) = r (1.1.2) A T A A A Din Ai TA, cum Ti sunt ineversabile găsim A T 1Ai şi deci: i i r i min (r i ,r ) = min (m,r i ) = r i (1.1.3) A T A A A Din (1.1.2) şi (1.1.3) rezultă că r i = r , i=1,3. A A Definiţia (1.1.3). Două matrice A şi B care se obţin una din alta prin transformări elementare se numesc echivalente ( în privinţa rangului) şi scriem A B. 1.1.3. FORMA GAUSS – JORDAN A UNEI MATRICE Fie A matricea unui sistem de m ecuaţii algebrice liniare cu n necunoscute (m < n). Definiţia 1.1.4. Matricea A se spune că are forma Gauss-Jordan dacă conţine r (r m) coloane ale matricei unitate de ordinul m. Dacă aceste coloane sunt j , j2,....., jr atunci forma matricei este 1 ....1.... 0 .... 0 a'k .... a'n  .... 0 ....1.... 0a2k .... a2n  1 1 ' '     ...................................... ' '       A  .... 0 .... 0 ....1ark .... arn (1.1.4) 0.. 0.... 0 .... 0 0 ....... 0 ......................................     0.. 0.... 0 .... 0 0 ....... 0  9 Teorema 1.1.2. Orice matrice nenulă poate fi adusă la forma Gauss-Jordan, printr-un număr finit de transformări elementare. Demonstraţie. Fie deci matricea A de tipul m x n a unui sistem liniar de m ecuaţii cu n necunoscute. 11 12 1j 1n a21 a22 1j a2n a a a a a ai1 ai2 aij ain am1 am2 amn amn Fie 11 0. Dacă 11 0, atunci printr-o transformare de tipul (T3) putem aduce în locul a a lui 11 un element i1 0; i 2,m care există deoarece prima coloană are cel puţin un element a a diferit de zero, în caz contrar sistemul nu ar avea n necunoscute. Printr-o transformare de tipul (T1) cu 1 realizăm în poziţia (1,1) un element egal 11 a cu unu. Prin m 1 transformări de tipul (T1) şi (T2) în care linia întâi este succesiv inmulţită cu numere convenabile: a21, a31, a1,..., am1 şi adunată de fiecare dată la liniile: a doua, a treia, i . . . a i-a linie, . . . a m-a linie, realizăm pe prima coloană toate elementele nule cu excepţia primului element care este unu. Dacă dorim să realizăm pe coloana j toate elementele nule cu excepţia elementului din poziţia (i,j) care să fie egal cu unu procedăm astfel: a) ne asigurăm că ij 0 a b) aplicăm o transformare de tipul (T1) cu 1 , obţinând în poziţia (ij) un element ij a egal cu unu. c) Aplicăm succesiv m 1 transformări de tipul (T1) şi (T2), înmulţinâd linia i cu: a j, a2 j , , amj şi adunând-o la liniile: întâia, a doua, . . . , a m-a linie obţinem 1 pe coloana j toate elementele nule cu excepţia elementului din poziţia (i,j) care devine unu. Algoritmul se opreşte când pe liniile r 1,m nu mai avem elemente nenule. Elementul nenul aij dat se va numi pivot, iar ansamblul transformărilor elementare 10 necesare transformării matricei A într-o matrice cu coloana i din matricea unitate de ordinul m în coloana j se va numi pivotajul cu element pivot aij . Din forma Gauss-Jordan a matricei A putem observa că rangul matricei este r, acesta fiind de ordinul minorului format cu primele r linii şi coloanele j , j2,...jr . Toţi minorii de 1 ordinul r1 sunt nuli. În practică nu este neapărat necesar să obţinem coloanele matricei unitate în ordinea lor naturală, numărul lor ne va da rangul. Putem alege pivot orice element diferit de zero oriunde s-ar găsi el de preferinţă cel care este egal cu 1. Exemplul 1.1.1. Să se aducă la forma Gauss-Jordan şi să se determine rangul matricei: ① 2 2 2 1 Luăm pivot elementul încercuit. Prima linie 1 2 1 1 2 înmulţită cu 1 se adună la a doua, apoi cu –2 şi o 2 4 4 3 1 adunăm la a treia obţinând matricea : 1 2 2 2 1 0 0 ③ 3 1 0 0 0 1 1 1 2 0 0 5/3 Luăm pivot elementul încercuit. Linia a doua o înmulţim cu 1/3 şi apoi înmulţită cu –2 o adunăm la prima obţinând matricea: Luăm pivot elementul încercuit. Linia a treia o 0 0 1 1 1/3 0 0 0 ① 1 înmulţim cu –1, apoi înmulţită cu –1 o adunăm la a doua şi obţinem: 1 2 0 0 5/3 0 0 1 0 2/3 0 0 0 1 1 Rangul matricei este egal cu 3, egal cu numărul maxim de coloane ale matricei unitate. 11 Exemplul 1.1.2. Să se determine rangul matricei: 2 1 1 3 1 2 1 2 1 1 2 1 Făcând pivotajele cu elementele pivot încercuite obţinem succesiv: 2 ① 1 3 2 1 1 3 1 2 1 2 3 0 1 4 1 1 2 1 3 0 ① 2 5 1 0 5 0 1 0 0 1 0 0 1 ⑥ 0 0 6 1 0 0 1 0 1 0 0 3 0 1 2 0 0 1 1 0 0 1 1 Rangul matricei considerate este 3. 1.1.4 INVERSA UNEI MATRICE Fie A o matrice pătrată de ordinul m şi B= A Im matricea obţinută prin alăturarea matricei unitate de ordinul m la matricea A. Se cunoaşte că se numeşte matrice inversă a unei matrice pătrate A, o matrice, notate cu A-1, astfel încât: AA-1=AA=Im (1.1.6) Presupunem pentru moment că există matricea A-1. Înmulţind la stânga matricea B cu A-1 obţinem: B=A 1B= A 1A A 1Im = Im A 1 De la matricea B putem ajunge la matricea B cu ajutorul transformărilor elementare. Făcând pivotaje în matricea A în vederea obţinerii celor m coloane ale matricei unitate de ordinul m, în dreapta barei vom găsi tocmai matricea A-1. Dacă coloanele nu sunt în ordinea lor naturală, putem schimba liniile între ele obţinând în stânga barei chiar matricea unitate. Dacă nu reuşim să obţinem coloane ale matricei unitate în număr egal cu ordinul matricei A, atunci matricea nu are inversă. 12 Exemplul 1.1.3. Să se găsească inversa matricei: 1 2 3 A= 1 1 2 1 1 1 considerând în matricea B pivotaje cu elemente pivot încercuite avem : ① 2 3 1 0 0 1 2 3 0 0 0 B 1 1 2 0 1 0 0 1 ① 1 1 0 1 1 1 0 0 1 0 1 2 1 0 1 1 1 0 ∣ 2 3 0 1 0 0 1 1 1 0 1 1 1 1 0 0 0 1 0 1 1 0 ① 0 1 2 1 0 1 0 1 2 1 1 0 0 1 1 1 0 1 0 1 2 1 0 0 1 0 1 1 Inversa matricei A este: 1 1 1 A 1= 1 2 1 0 1 1 Exemplul 1.1.4. Să se găsească inversa matricei: 1 1 2 A= 2 2 1 3 1 1 Vom avea: ① 1 2 1 0 0 1 1 2 1 0 0 B 2 2 1 0 1 0 0 ④ 5 2 1 0 3 1 1 0 0 1 0 4 5 3 0 1 1 0 3/4 1/2 1/4 0 0 1 5/4 1/2 1/4 0 0 0 0 1 1 1 Matricea nu admite inversă. 13 1.2. SISTEME DE ECUAŢII LINIARE Un sistem de ecuaţii algebrice liniare cu m ecuaţii şi n necunoscute (m < n) se prezintă sub forma: 11 1 12x2 1nxn b , a x a a 1 a x a x a x b , (1.2.1) . . . . . . . . . . . . . . . . . . . . . . . . . 21 1 22 2 2n n 2 x am1 1 am2x2 amnxn bm, unde aj,b (i=1,m ; Dacă există i i i b 1,n) sunt numere reale. 0 sistemul (1.2.1) este neomogen, în caz contrar (b 0,( )i 1,m) i sistemul se numeşte omogen. Fie B matricea lărgită a sistemului (1.2.1), de forma: 11 12 1n b a a a 1 a21 a22 a2n b2 (1.2.2) am1 am2 amn bm Să presupunem că rangul matricei A a sistemului este r (r m) . Prin transformări elementare putem aduce matricea A la forma Gauss-Jordan. Obţinem astfel matricea B , de forma: 1 ... 0 0 1jr 1 1n b 0 ... 1 0 a2 jr 1 a2n b2 . ... . ... . . ... . . a a 1 B 0 ... 0 ... 1 arjr 1 arn br (1.2.3) 0 ... 0 ... 0 0 ... 0 r 1 b . ... . ... . . ... . . 0 ... 0 ... 0 0 ... 0 bm Sistemul corespunzător matricei B este: 14 a a b 1 r xj xj 1 1 1jr 1 xjr 1 2 jr 1 jr 1 a x 1nxn 1 a2nxn b2 xjr a jr 1 xjr 1 arnxn br (1.2.4) 0 b 1 0 bm Deoarece matricele B şi B sunt echivalente, sistemele (1.2.1) şi (1.2.4) vor fi echivalente (adică vor avea aceleaşi soluţii). Este suficient pentru aceasta să observăm că transformările elementare aplicate liniilor matricei B pentru a fi adusă la forma B au drept consecinţe asupra sistemului de ecuaţii: (T1) – înmulţirea unei ecuaţii cu un număr diferit de zero; (T2) – adunarea la o ecuaţie a unei alte ecuaţii (eventual înmulţită cu un scalar nenul); (T3) – schimbarea a două ecuaţii între ele. Evident, aceste operaţii aplicate unui sistem de ecuaţii îl transformă în unul echivalent. Teorema 1.2.1. Condiţia necesară şi suficientă ca sistemul (1.2.1) să fie compatibil este ca b 0, i r 1,m. i Demonstraţie. Să arătăm mai întâi necesitatea condiţiei. Presupunem că sistemul (1.2.1) este compatibil. Aceasta înseamnă că sistemul (1.2.4) este compatibil. Din compatibilitatea acestuia din urmă rezultă i 0 pentru i r 1,m. b Să arătăm acum suficienţa condiţiei. Presupunem că b 0, i r 1,m. Rezultă atunci că i sistemul (1.2.4) este compatibil admiţând n r soluţii şi deci sistemul (1.2.1) este compatibil. În baza acestei teoreme uşor pot fi demonstrate teoremele lui Kronnecker-Capelli şi Rouché. Pentru r m n se obţin rezultatele cunoscute de la sistemele de n ecuaţii cu n necunoscute. 15 Exemplul 1.2.1. Să se atudieze compatibilitatea sistemului şi în caz afirmativ să se găsească soluţia lui. 1 x2 x3 x4 1 x x2 2x3 2x4 2 x 1 3 x x2 3x3 3x4 5 x 3x2 5x3 5x4 3 1 1 Matricea lărgită a sistemului, după o serie de pivotaje cu elemente pivot încercuite, devine: B ① 1 1 1 1 1 1 2 2 2 3 1 3 3 5 1 3 5 5 3 1 1 1 1 1 1 0 1/2 1/2 3/2 0 ② 3 3 1 0 1 3/2 3/2 1/2 0 4 6 6 2 0 0 0 0 0 0 4 6 6 2 0 0 0 0 0 Sistemul considerat este echivalent cu sistemul: 1 1/2x3 1/2x4 3/2 x2 3/2x3 3/2x4 1/2 x . care este compatibil având soluţia:  x 3/21/2 1/2,  x2  1/23/2 3/2, x3 , x4  , ,  1   Exemplul 1.2.2. Să se studieze compatibilitatea sistemului: 1 2x2 3x3 4x4 2 x x2 x3 x4 3 x 1 2x 5x2 6x3 x4 10 1 16 Matricea lărgită a sistemului după o serie de pivotaje cu elemente pivot încercuite devine: ① 2 3 4 2 1 2 3 4 2 B 1 1 1 1 3 0 ③ 4 3 5 3 5 6 1 10 0 9 12 9 14 1 0 1/3 2 4/3 0 1 4/3 1 5/3 0 0 0 0 1 Cum b 1 0, sistemul este incompatibil. 3 1.2.1. EXPLICITAREA SISTEMELOR DE ECUAŢII LINIARE Fie sistemul (1.2.1). Presupunem că rangul matricei A este m < n. Definiţia (1.2.1). Vom spune că sistemul (1.2.1) este explicitat în raport cu un grup de m variabile (necunoscute), dacă în matricea A a sistemului coloanele acestor variabile sunt cele m coloane ale matricei unitate de ordinul m. Deoarece cu cele n variabile ale sistemului putem forma Cm grupuri diferite de câte m n variabile, rezultă că un sistem liniar, care poate fi explicitat în raport cu cel puţin un grup de m variabile, va avea cel mult Cm forme explicite. Acest număr maxim va fi atins dacă explicitarea n poate avea loc în raport cu oricare dintre cele Cm grupuri diferite de câte n variabile şi toate n formele explicite sunt distincte. Într-un sistem liniar explicitat sunt două tipuri de variabile (necunoscute): unele, ai căror coeficienţi formează coloanele matricei unitate, altele ai căror coeficienţi formează celelalte coloane ale matricei sistemului. Ca de obicei, primele vor fi denumite variabile principale, iar celelalte variabile secundare. În teoria programării liniare variabilele principale mai sunt denumite bazice, iar variabilele secundare nebazice. Definiţia (1.2.2). Numim soluţie de bază a unui sistem de ecuaţii liniare, o soluţie particulară obţinută dintr-o formă explicită, prin egalarea cu zero a variabilelor secundare. 17 i a i a m Din definiţia dată rezultă că dacă sistemul este compatibil nedeterminat, atunci are cel puţin o soluţie de bază şi cel mult Cm soluţii de bază. Evident, sistemele incompatibile şi cele n compatibil determinate nu au soluţii de bază. Dacă o soluţie de bază are exact m componente nenule, atunci soluţia de bază se numeşte nedegenerată, iar dacă are mai puţin de m componente nenule se numeşte soluţie de bază degenerată. Soluţiile de bază ale unui sistem liniar se clasifică şi după un alt criteriu. Dacă toate componentele unei soluţii de bază au valori nenegative ( 0), soluţia de bază se numeşte admisibilă. În caz contrar, soluţia de bază se numeşte neadmisibilă (are componente strict negative). Să considerăm că r = m. În acest caz din (1.2.4) obţinem o soluţie de bază: xj b , xj2 b , . . . xji b ,..... xjm bm , restul până la n egale cu zero. Dorim ca variabila 1 1 2 i principală xji să devină secundară iar variabila secundară xj ( j j , . . . ,jr ) să devină 1 principală. Pentru aceasta fie şi ij 0. Efectuând în matricea B un pivotaj cu element pivot a aij găsim o altă soluţie de bază: xj b a j i ,xj ij 1 1 b 1 a 2 b2 a2 j b ,. . .,xj ij b , . . .xj bm ij i mj b a a ij şi restul până la n egale cu zero. Putem scrie variabilele principale concentrat sub forma: b xjk xj bk b i aij i kj a a ij , pentru k j (1.2.5) , pentru k j Exemplul 1.2.3. Să se determine toate formele explicite şi soluţiile de bază corespunzătoare pentru sistemul: 1 x2 x3 0 1 2x2 x3 x4 2 x x Matricea lărgită a sistemului este: 18 1 x x 1 x x 1 x x 2 1 B 1 1 1 0 0 1 2 1 1 2 Pivotajele: a) ① 1 1 0 0 1 1 1 0 0 1 0 1 0 0 1 2 1 1 2 0 1 0 ① 2 0 1 0 1 2 b) ① 1 1 0 0 1 1 1 0 0 1 0 1 1 2 1 2 1 1 2 0 ① 0 1 2 0 1 0 1 2 c) 1 ① 1 0 0 1 1 1 0 0 0 1 0 1 2 1 2 1 1 2 1 0 ① 1 2 1 0 1 1 2 d) 1 1 ① 0 0 1 1 1 0 0 1 1 1 0 0 1 2 1 1 2 0 1 0 ① 2 0 1 0 1 2 e) 1 ① 1 0 0 1 1 1 0 0 1 1 1 0 0 1 2 1 1 2 1 0 1 ① 2 1 0 1 1 2 conduc la formele explicite în raport cu grupurile: (x ,x2),(x ,x4),(x2,x3),(x2,x4),(x3,x4), 1 1 a) x + x3 x4 2 4 2 2 b) x + x3 2 2 x4 2 c) x2 1 + x 2 x3+ x4 2 4 d) x +x2 x3 1 3 0 x4 2 f) 1 1 x x x 2x2 x3 4 x 0 2 Soluţiile de bază vor fi: (S ) x 2, x2 2, x3 0, x4 0; soluţie de bază nedegenerată admisibilă 1 1 (S2) x 2, x2 0, x3 0, x4 2; soluţie de bază nedegenerată admisibilă 1 (S3) x 0, x2 (S4) x 0, x2 1 (S5) x 0, x2 1 2, x3 2, x4 0, x3 0, x4 0, x3 0, x4 0; soluţie de bază nedegenerată admisibilă 2; soluţie de bază degenerată admisibilă 2; soluţie de bază degenerată admisibilă Sistemul nu poate fi explicitat în raport cu grupul de variabile (x ,x3) şi admite doar 1 patru soluţii de bază distincte din numărul maxim de C2  6 posibilităţi. 4 19 US2 Spații liniare Noțiunea de vector și cea de spațiu liniar sunt noțiuni centrale în algebra liniară. Utilizarea acestora pentru rezolvarea problemelor de optimizare este cât se poate de naturală. În prezenta unitate de studiu se prezintă noțiunea de spațiu liniar cu exemplul cel mai potrivit (adică spațiul n ) dar și cu alte exemple. De asemenea se aplică elemente din US1 (Transformări elementare) pentru a stabili relația de dependență (independență) liniară a unor familii de vectori. De asemenea se utilizează sistemele liniare (abordate și în US1) pentru a stabili concret relațiile de dependență liniară. Se introduce în continuare noțiunea de bază de vectori, precum și de coordonate ale unui vector într-o bază. Schimbarea de coordonate, atunci când se schimbă vectorii bazei, este următorul subiect abordat în material. Cei care și-au însușit complet US1 nu vor avea dificultăți pentru parcurgerea acestei unități de studiu. Mai întâi trebuie însușite corect noțiunile de bază (vector, bază, coordonate). Ulterior trebuie aplicate teoremele de caracterizare a dependenței și independenței (folosind efectiv transformările elementare). Pentru determinarea coordonatelor unui vector se vor rezolva (cu transformări elementare) sisteme de ecuații compatibile determinate. Însușirea noțiunilor teoretice presupune un studiu de patru ore. Aplicații pentru dependență și independență solicită încă două ore. Deteminarea coordonatelor (rezolvarea sistemelor liniare) cu aplicațiile corespunzătoare mai adaugă încă două ore la timpul de studiu. În capiutolul acesta se prezintă un algoritm de calcul a noilor coordonate ale unui vector când unul dintre vectorii bazei este înlocuit. Însușirea lui și aplicarea cu ușurință înseamnă încă două ore de studiu. În total US2 necesită zece ore de studiu. Pentru facilitarea studiului sunt utile din bibliografie pozițiile [5] și [6]. * * * Înzestrând o mulţime oarecare cu o serie de operaţii având anumite proprietăţi obţinem noţinea de structură algebrică. Unele tipuri de structuri algebrice au fost studiate în liceu (monoid, grup, inel, corp). Spaţiul liniar se defineşte în acelaşi mod, elementele sale (vectorii) generalizând noţiunea de vector liber utilizat în fizică sau matematică. 20 În orice spaţiu liniar se introduce noţiunea de independenţă liniară precum şi cea de bază. Vom introduce aceleaşi noţiuni dar numai în cazul particular al spaţiului liniar n (care apare cel mai des în contextul aplicaţiilor economice). În studiul problemelor economice se impune problema determinării valorii maxime sau a valorii minime a unor funcţii. Un ajutor preţios în acest scop îl constituie formele pătratice care generalizează funcţia de gradul al doilea. 21 2.1. DEFINIŢIA SPAŢIULUI LINIAR Fie A o mulţime oarecare şi K un corp având elementele neutre 0K respectiv 1K. Presupunem că pe A se pot defini două operaţii astfel: a) ( )x,y A, ( ) un element notat cu x y A astfel încât (x,y) x y. b) ( ) x A, ( ) K, ( ) un element notat x A, astfel încât ( ,x) x. Prima operaţie, de tip aditiv, este o operaţie internă iar cea de doua, de tip multiplicativ, este o operaţie externă pentru A. Definiţie. Mulţimea A formează un spaţiu liniar peste corpul K dacă sunt satisfăcute următoarele proprietăţi: L1) ( )x, y, z A, x (y z) (x y) z, L2) ( ) x, y A, x y y x, L3) ( )0A A astfel încât x 0A x, ( ) x A, L4) ( ) x A, ( ) x A astfel încât x ( x) 0A , L5) ( ) K, ( ) x,y A, (x,y) x y, L6) ( ) K, ( ) K, ( ) x A, ( ) x x y , L7) ( ) , K, ( )x A, ( x) ( x) ( )x, L8) ( )x A, 1K x x. Elementele lui K se numesc scalari, iar cele ale lui A, vectori. Prima operaţie se numeşte adunarea vectorilor iar cea de a doua înmulţirea vectorilor cu scalari din K. Atunci când nu este pericol de confuzie vectorul nul 0A se notează simplu 0. De asemenea, elementele neutre ale lui K se scriu 0 şi 1 (în imensa majoritate a aplicaţiilor K= ). Consecinţe. Din axiome decurg următoarele proprietăţi: a) ( ) x A, 0x 0, Într-adevăr, din L6) obţinem ( 0) x x 0x, adică x x 0x. Axioma L3 ) conduce la 0x 0. b) ( ) K, 0 0. Plecăm de la L5) : (x 0) x 0. Obţinem x x 0. 22 Aceeaşi axiomă L3) arată că: 0 0. c) ( ) x A, ( 1) x x Folosim a), adică 0x 0 şi înlocuim 0 prin 1 ( 1). Din L6) obţinem: 1x ( 1)x 0. L8) conduce la x ( 1)x 0, iar L4) arată că ( 1)x 0. Exemplu: Considerăm A n x x x (x ,x2, ,xn)/ x şi K= . 1 i Adunarea elementelor din A revine la adunarea matricelor cu o singură coloană sau cu o singură linie. ( )x (x ,x2, ,xn), y (y , y2, , yn) din n , 1 1 x y (x y , x2 y2, xn yn) n 1 1 În mod similar produsul unui scalar din cu un element al lui n se defineşte prin: x (x ,x2, ,xn) ( x , x2,......, xn). 1 1 2.2. DEPENDENŢĂ ŞI INDEPENDENŢĂ LINIARĂ Fie A un spaţiu liniar peste corpul K şi fie xi A, i K, i 1,h. Definiţie 2.2.1. Se numeşte combinaţie liniară a vectorilor xi cu scalarii i expresia : E 1x1 2x2 hxh . (2.2.1) Consecinţa a) din definiţia spaţiului liniar ne arată că atunci când i 0,( ) i 1,h E 0. Reciproca acesteia nu este întotdeauna adevărată. Definiţie 2.2.2. Vectorii xi A, i 1,h se numesc liniar independenţi dacă: 1x1+ 2x2+ + hxh 0 1 2 h 0. În caz contrar, dacă E=0 cu măcar un scalar nenul, vectorii xi se numesc liniari dependenţi, adică: 23 1x1+ 2x2+ + hxh 0 1 2 h 0 Exemplu: Fie spaţiul liniar A n n,1( ) şi vectorii: x1 (x 1,x 2, ,x n)t,x2 (x21,x22, ,x2n)t, ,xh (xh1,xh2, xhn)t vectori din A. 1 1 1 Condiţia E 0 conduce la sistemul liniar şi omogen: 1 11 2x21 hxh1 0 1 12 2 22 h h2 x x x x 0 . . . . . . . . . . . . . . . . . . . . . . . (2.2.3) 1 1n 2x2n hxhn 0 x Vectorii sunt liniari independenţi dacă sistemul (2.2.3) admite numai soluţia banală şi liniar dependenţi în caz contrar. Din teoria sistemelor omogene independenţa liniară este caracterizată prin: rang M r = h (2.2.4) M unde matricea M n,h ( ) este formată prin componentele vectorilor xi . Dacă (2.2.3) are soluţii nebanale, ele se exprimă prin h r parametri reali. Înlocuite în M (2.2.2) aceste soluţii conduc la una sau mai multe relaţii de dependenţă liniară. Pe baza acestui exemplu vom realiza un mare număr de aplicaţii. Proprietăţile familiilor de vectori liniar independenţi şi al familiilor de vectori liniar dependenţi sunt date în următoarele teoreme: Teorema 2.2.1. O mulţime de vectori liniar independenţi nu poate conţine vectorul nul. Demonstraţie: Presupunem prin absurd că mulţimea de vectori independenţi x1,x2, xi, ,xh ar conţine vectorul nul şi că xi 0. Combinaţia liniară : E 0 x1 0 x2 1 xi hxh 0 cu scalarul i nenul. Contradicţia justifică teorema. Teorema 2.2.2. Orice submulţime a unei mulţimi de vectori liniar independenţi este formată tot din vectori liniar independenţi. 24 Demonstraţie: Tot prin reducere la absurd presupunem că deşi x1,x2, ,xh sunt liniar independenţi, submulţimea x1 ,xi2 , ,xik , k h este formată din vectori liniar dependenţi. i Atunci, egalitatea: i xi + i2 xi2 ik xik 0 1 1 are loc cu măcar un scalar nenul. Rezultă că: i1xi1 + i2 xi2 ik xik 0 x 0 j 1,h j j I cu măcar un scalar nenul ceea ce contrazice ipoteza teoremei. Teorema 2.2.3. Într-o mulţime de vectori liniar dependenţi cel puţin un vector se scrie ca o combinaţie liniară a celorlalţi. Demonstraţie: Definiţia vectorilor liniar dependenţi arată că: 1x1+ 2x2+ ixi + hxh 0 cu măcar un scalar nenul. Presupunem că acesta este i . Cum orice spaţiu liniar este grup abelian, regulile de calcul dintr-o astfel de structură permite să scriem: ixi 1x1 2x2 .......... ixi 1 i 1xi 1 ........... hxh În corpul K, i este inversabil iar dacă i 1 este inversul său xi 1 i 1 1 x ...... 1 i i 1xi 1 1 i i 1xi 1 ...... 1 i hxh (2.2.5) şi teorema este demonstrată. Când corpul de scalari este inversul i 1 al lui i este 1/ i şi obţinem: xi 1 x1 2 x2 i i i 1 i 1 i 1 i 1 i i x x h xh (2.2.6) i Teoremele prezentate au consecinţe importante cu valoare teoretică şi practică. Astfel, modificarea componenţei unei mulţimi de vectori liniar independenţi prin suprimarea unora conduc totdeauna la o familie de vectori liniar independenţi. Dacă adăugăm noi vectori unei astfel de familii, rezultatul nu este previzibil, noua mulţime poate fi dependentă sau independentă. Din (2.2.3) rezultă că (2.2.4) nu poate fi realizată dacă h>n, deci numărul maxim de vectori liniar independenţi din n este n. 25 Existăefectivmulţimi devectoriliniarindependenţi. Unadintreacesteaeste: F= e1,e2,…,en (2.2.7) unde: e1 (1,0,0, ,0)t, e2 (0,1,0, ,0)t, ,en (0,0,0, ,1)t . Dacă prin transformări elementare modificăm matricea M a vectorilor mulţimiiF , obţinem alte mulţimi formate tot din n vectori liniar independenţi. Numărul acestor mulţimi este infinit. 2.3. BAZE DE VECTORI. COORDONATELE UNUI VECTOR ÎNTR-O BAZĂ Am luat în studiu în paragraful precedent familiile maximale de vectori liniar independenţi dintr-un spaţiu liniar. Rolul acestora va fi precizat în cele ce urmează. Definiţie. Familia B u1,u2, ,uk formează o bază în spaţiul vectorial V dacă: a) u1,u2, ,uk sunt liniar independenţi; b) ( ) u V, mulţimea F u,u1,u2, ,uk , este formată din vectori liniar dependenţi. Mulţimea B este aşadar o mulţime maximală de vectori liniar independenţi. Când spaţiul V este n , această mulţime conţine n vectori, deci k = n. Spaţiul n are cel puţin o bază: cea formată din e1,e2,…,en numită bază canonică. Aşadar: Bc e1,e2,…,en Din aceasta prin procedeul indicat anterior putem obţine o infinitate de baze toate având n vectori. Numărul n, al vectorilor oricărei baze din n se numeşte dimensiunea acestui spaţiu. Condiţia b) din definiţia bazei şi una din proprietăţile fundamentale ale unei familii de vectori liniar dependenţi ne permit să formulăm şi să demonstrăm: 26 Teorema 2.3.1. Orice vector din n se exprimă în mod unic printr-o combinaţie liniară a vectorilor unei baze B. Demonstraţie. Teorema 2.2.3 ne arată că cel puţin un vector al mulţimii u,u1,u2, ,un în care u1,u2, ,un formează bază în n se exprimă printr-o combinaţie liniară a celorlalţi. Acest lucru este valabil şi pentru u, deoarece în caz contrar, din: u 1u1 2u2 nun 0 (2.3.1) în care 0 ar rezulta: 1u1 2u2 nun 0 (2.3.2) şi unul din scalari fiind nenul s-ar contrazice condiţia a) din definiţia bazei. Aşadar 0 şi obţinem: u 1 u1 2 u2 n un (2.3.3) Notăm: 1 1 , 2 2 , , n n şi scriem: u 1u1 2u2 nun (2.3.4) Prin reducere la absurd demonstrăm că scalarii i din (2.3.4) sunt unic determinaţi. În caz contrar, are loc şi relaţia: u 1u1 2u2 nun (2.3.5) Din (2.3.4)-(2.3.5), obţinem: ( 1 1)u1 ( 2 2)u2 ( n n)un 0 (2.3.6) Cum u1,u2, ,un sunt liniar independenţi, toţi scalarii din (2.3.6) sunt nuli, deci: 1 1, 2 2, , n n ceea ce este în contradicţie cu ipoteza. Teorema este demonstrată. 27 3 Observaţie. Scalarii i, (i 1,n) se numesc în mod uzual coordonatele vectorului u în baza B u1,u2, ,un . Unicitatea lor rezultă şi din teoremă dar şi din faptul că (2.3.4) conduce la un sistem Cramer (cu soluţie unică). Vom nota mulţimea coordonatelor vectorului u în baza B cu uB , adică: uB ( 1, 2, , n )t (2.3.7) Un calcul elementar ne arată că în baza canonică Bc e1,e2,…,en coordonatele oricărui vector u coincid cu componentele sale, adică: u (u1,u2, ,un)t u1e1 u2e2 unen Exemplu: Să se arate că vectorii u1 (1,1,1)t, u2 (1,2,1)t, u3 (1,1,3)t formează o bază în 3 şi să se determine coordonatele lui u (6,1,8)t în această bază. Relaţia (2.3.4) devine: u 1u1 2u2 3u3 şi conduce la sistemul: 1 1 1 2 2 2 2 3 3 3 6 1 8 Prin transformări elementare de linii în matricea lărgită a sistemului obţinem: 1 1 1 6 1 1 1 6 1 2 1 1 0 1 0 5 1 1 3 8 0 0 2 2 1 0 1 11 1 0 0 10 0 1 0 5 0 1 0 5 0 0 2 2 0 0 1 1 care corespunde sistemului cu r(M) 3 ) deci (u1, u2, u3 formează baza). 1 2 3 10 5 1 acestea fiind valorile coordonatelor căutate. 28 2.4. SCHIMBAREA COORDONATELOR UNUI VECTOR LA SCHIMBAREA BAZEI Dacă în locul bazei B u1,u2, ,un a spaţiului liniar n considerăm o altă bază B v1,v2, ,vn având acelaşi număr de vectori, n, ne aşteptăm ca în relaţia (2.3.4) coordonatele să se schimbe. Dorim să aflăm legea de schimbare a acestora în cazul general şi într-un caz particular. Să considerăm că în baza B (numită baza veche) vectorul u are coordonatele 1, 2, , n , deci: u 1u1 2u2 nun (2.4.1) şi că acelaşi vector are în baza B (numită baza nouă) coordonatele 1, 2, n , adică: u 1v1 2v2 nvn (2.4.2) Pe de altă parte fiecare din vectorii vi ai bazei noi are în baza veche coordonatele aji , unic determinate vi a1iu1 a2iu2 aniun jn1ajiuj (2.4.3) Înlocuim (2.4.3) în (2.4.2) şi comparăm cu (2.4.1). Obţinem: i ai1 1 ai2 2 ain n, i 1,n (2.4.4) Dacă notăm cu A matricea coordonatelor vectorilor din baza nouă faţă de vectorii din baza veche (numită matrice de trecere): a11 a21 an1 A a12 a22 an2 a1n a2n ann cu uB =( 1, 2, , n )t şi cu uB' =( 1, 2, , n )t atunci (2.2.4) se scrie matricial : uB AtuB (2.4.5) 29 2 Cum vectorii lui B sunt liniar independenţi, A este inversabilă şi din (2.4.5) putem obţine relaţia: uB At 1 uB (2.4.6 ) Relaţiile (2.4.5) şi (2.4.6) sunt formulele de schimbare a coordonatelor la schimbarea bazei. Exemplu: Se consideră vectorii u1 (1,2,1)t, u2 (2,1,0)t, u3 (1,1,0)t din 3 : a) Să se arate că formează o bază în 3 b) Să se determine coordonatele lui u (1,1, 2)t în această bază. c) Să se determine coordonatele lui u în baza formată din v1 u1 u2, v2 u1 u2 u3, v3 u2 2u3 Soluţie: a) Matricea formată cu componentele vectorilor u1,u2,u3 are rangul trei, egal cu dimensiunea lui 3 , aşadar cei trei vectori formează o bază. b) Vectorul u se scrie u 1u1 2u2 3u3 care conduce la sistemul 1 2 2 3 1 2 3 1 2 1 1 cu soluţia: 1 2, 2 2, 3 7, deci : uB ( 2, 2,7)t,B u1,u2,u3 c) Matricea de trecere de la B la B v1,v2,v3 este : 1 1 0 A 1 1 1 0 1 2 30 iar inversa lui At are forma: 3/5 2/5 1/5 (At ) 1 2/5 2/5 1/5 1/5 1/5 2/5 Folosind relaţia (2.4.6), obţinem în final: uB ( 1, 3,2)t . Observaţie. Dacă vectorii noii baze nu sunt exprimaţi de la început prin combinaţii ale vectorilor vechii baze, coordonatele lor (deci elementele matricei de trecere A) trebuie calculate separat. În această situaţie este mai simplu să calculăm direct uB folosind definiţia. În unele cazuri bazele B şi B diferă printr-un singur vector. Formula de calcul a noilor coordonate va avea evident o formă mai simplă. Să presupunem că în baza B, vectorul ui s-a înlocuit printr-un vector vj care are coordonatele a1j,a2j, ,aij, ,anj , adică: vj a1ju1 a2ju2 aijui anjun (2.4.7) Ceilalţi vectori ai vechii baze rămânând neschimbaţi, matricea de trecere A devine: 1 0 0 0 0 1 0 0 a1j a2j aij anj 0 0 0 1 După determinarea lui (At ) 1 şi aplicarea formulei (2.4.6) deducem noile coordonate: i i , k k i akj, (k i) (2.4.8) ij ij a a Formulele (2.4.7), greu de memorat pot fi utilizate prin intermediul unui algoritm de calcul bazat pe transformări elementare ale liniilor matricei coordonatelor vectorilor u,u1,u2, ,ui, un şi vj . Ilustrăm mai jos modul de organizare a algoritmului. Etapa iniţială are forma : 31 a Baza u1 u2 … ui … un vj uB u1 1 0 … 0 … 0 a1j 1 u2 0 1 … 0 … 0 a2j 2 … … … … … … … … ui 0 0 … 1 … 0 aij i … … … … … … … … un 0 0 … 0 … 1 anj n Coloanele tabelului conţin coordonatele fiecărui vector în baza de pe prima coloană. Ele se presupun cunoscute. Ultima coloană reprezintă coordonatele lui u în baza veche B. Etapa finală are forma: Baza u1 u2 … vj … un vj uB u1 1 0 … u2 0 1 … a1j aij a2j aij … 0 0 1 1 … 0 0 2 2 i a 1j a ij i a 2j a ij … … … … … … … … vj 0 0 … 1 aij … 0 1 i i aij … … … … … … … … un 0 0 … anj aij … 1 0 n n i nj a ij şi conţine coordonatele tuturor vectorilor în noua bază. Se observă că linia i din ultima fază se obţine prin împărţirea liniei i din prima fază prin aij (care este nenul, altfel matricea de trecere A nu este nesingulară). 32 O linie oarecare k a matricei finale se obţine înmulţind linia i cu aij după care se adună cu linia k a matricei iniţiale. Ca efect al acestor operaţiuni, pe ultima coloană se obţin coordonatele i, k(k i), pe coloana vj se obţine coloana "i" a matricei unitate In iar coloana lui vj se modifică în mod corespunzător. Algoritmul poate fi aplicat înlocuind cauzele cu efectele şi invers. Astfel, putem să aplicăm transformările elementare de linii în matricea fazei iniţiale pentru a transforma coloana lui vj în coloana vectorului ui din faza iniţială. Unii numesc această operaţie “Lema substituţiei” şi aplică diverse metode mnemotehnice de utilizare. Noi o considerăm o simplă aplicare a transformărilor elementare şi o folosim succesiv înlocuind toţi vectorii bazei iniţiale B prin vectorii bazei B . Exemplu. În baza canonică Bc e1,e2,e3 se consideră vectorul u (1,1, 1)t . Să se determine coordonatele sale în baza B v1,v2,e3 unde v1 (1,2,3)t,v2 (1, 1,1)t,v3 (0,0,1)t . Baza e1 e2 e3 e1 1 0 0 e2 0 1 0 e3 0 0 1 v1 1 0 0 e2 2 1 0 e3 –3 0 1 v1 1/3 1/3 0 v2 2/3 –1/3 0 e3 –5/3 –2/3 1 v1 v2 u 1 1 1 2 –1 1 3 1 –1 1 1 1 0 –3 –1 0 –2 –4 1 0 2/3 0 1 1/3 0 0 –10/3 Operaţiile – extrem de simple – au urmărit crearea matricei unitate pe coloanele vectorilor noii baze B v1,v2,e3 . 33 În prima etapă am completat tabelul cu coordonatele – coincizând cu componentele – tuturor vectorilor în baza canonică Bc. S-au aplicat transformări elementare de linii pentru a obţine pe coloana lui v1 prima coloană a matricei unitate I3 . În ultima etapă, prin transformări elementare de linii se obţin pe coloana lui v2 cea de a doua coloană a matricei unitate I3 . Coordonatele lui u în noua bază sunt înscrise pe ultima coloană şi deci: u 2 v1 1 v2 10e3 . 3 3 3 Observaţie. Chiar dacă pare prea complicat în comparaţie cu aplicarea definiţiei în determinarea coordonatelor lui u, algoritmul are avantajul principal că în fiecare etapă coloanele conţin coordonatele în baza curentă. În plus, ca avantaj secundar este posibilitatea de a-l folosi în probleme de programare liniară. El va fi aplicat consecvent în capitolele care urmează. 34 US3 Forme pătratice Formele pătratice constituie un caz particular de funcții definite pe un spațiu liniar. Ele sunt întâlnite în numeroase capitole de matematici economice (chiar unele modele de creștere economică utilizează formele pătratice). Dedicăm această unitate de studiu formelor pătratice făcând mai întâi o prezentare a noțiunii, după care facem o clasificare (după semn) a acestora. Caracterizarea formei se obține cu ușurință după aducerea la expresia canonică prin una dintre metodele devenite clasice (metoda lui Jacobi și metoda lui Gauss). Sunt date aplicații în acest sens, accentuând aplicarea acelorasi transformări elementare de linii în matricea asociată formei. Studentul trebuie să-și însușească mai întâi noțiunile fundamentale. Ulterior trebuie să învețe să asocieze unei forme pătratice o matrice simetrică și reciproc. Ultimul lucru este să folosească metoda transformărilor elementare pentru a aplica metoda lui Gauss de obținere a expresiei canonice (sau a metodei lui Jacobi pentru a ajunge la același rezultat). Sunt necesare aproximativ șase ore pentru a finaliza studiul și a fixa, prin aplicații, cunoștințele (două ore pentru studiul teoretic, o ora pentru matricea asociată, două ore pentru metoda lui Gauss și o oră pentru metoda lui Jacobi). Sursele bibliografice sunt cele de pe pozițiile [3], [4], [5], [6] și [8]. * * * Considerăm în continuare spaţiul liniar n definit peste precum şi funcţia Q: n x n definită prin: Q(x) in1 jn1aijxixj, cu aij aji ( ) i, j 1,n. (3.1) Această funcţie se numeşte pătratică şi generalizează funcţiile de gradul al doilea omogene într-o singură variabilă. Matricea A n( ) care are elementele aij, i, j 1,n se numeşte matricea formei pătratice. Prin verificare directă se poate arăta că are loc formula: Q(x) xtAx. (3.2) 35 x x Exemplu: Să se scrie matricea corespunzătoare formei pătratice Q: 3 definită prin: Q(x) 2x2 x2 3x2 x1x2 2x2x3 . 1 2 3 Soluţie. Folosim (3.1) sau (3.2) elementul aij al matricei A fiind coeficientul lui xixj în expresia acesteia. Cum xixj xjxi şi aij aji în A coeficienţii produselor de variabile diferite se înjumătăţesc. Aşadar: 2 1/2 0 A 1/2 1 1 0 1 3 şi formula (3.2) se poate imediat verifica. Observaţie. Matricea A este pătratică şi simetrică conform condiţiei aij aji (nu este obligatoriu inversabilă). Dacă în matricea A numai elementele diagonalei principale pot fi nenule spunem că forma pătratică are expresia canonică. Expresia canonică a unei forme pătratice este: Q(x) 2 2 1 1 2 2 nxn , i . (3.3) 2 O formă pătratică Q este nedegenerată dacă expresia canonică are exact n termeni. Când numărul lor este mai mic forma se numeşte degenerată. Expresia canonică a unei forme pătratice are o importanţă deosebită în studiul semnului acesteia având în vedere că Q(x) . Există forme pătratice care păstrează acelaşi semn pentru orice argument x n . Nu este greu de verificat că: Q(0) 0. Următoarele definiţii pot fi relevante: Definiţie. Forma pătratică Q: n se numeşte pozitiv (negativ) definită dacă: Q(x) 0 (Q(x) 0) ( )x n, x 0. (3.4) 36 Definiţie. Forma pătratică Q: n se numeşte pozitiv (negativ) semidefinită dacă: Q(x) 0 (Q(x) 0) ( )x n (3.5) Orice formă pătratică definită este o formă pătratică semidefinită. Reciproca nu este adevărată. În cazul în care: ()x, y n a.i. Q(x) 0 si Q(y) 0 atunci forma pătratică se numeşte nedefinită ca semn. Dacă o formă pătratică este adusă la expresia canonică caracterul de formă pătratică pozitiv (negativ) definită (semidefinită) poate fi mult mai simplu apreciat. Astfel, dacă forma pătratică Q(x) definită de (3.1) nu este sub forma (3.3), putem obţine expresia canonică a acesteia: Q(y) 2 2 1 1 2 2 y y nyn , i (3.6) 2 unde yi , i 1,n sunt date de relaţiile liniare: 1 11 1 12x2 1nxn 2 21 1 22 2 2n y b x b b y b x b x b x . . . . . . . . . . . . . . . . . . . . . . . n ; cu bj (3.7) yn bn1 1 bn2x2 bnnxn x i Teoremă: Fie forma pătratică Q(x) definită de (3.1) a cărei expresie canonică asociată este dată de (3.6). Dacă:  a) i 0 , i 1,n b) i 0 , i 1,n  c) i 0 , i 1,n  d) i 0 , i 1,n  atunci Q(x) este pozitiv definită; si ( ) j 0atunci Q(x) este pozitiv semidefinită;  atunci Q(x) este negativ definită; si ( ) j 0atunci Q(x) este negativ semidefinită;  e) () i 0 si ()j 0 atunci Q(x) este nedefinită ca semn.  Vom da în continuare două metode de obţinere a expresiei canonice (3.6). a) Metoda lui Gauss. Se aplică transformări elementare asupra liniilor matricii A corespunzătoare formei pătratice, pentru a o aduce la forma (numită superior triunghiulară) de mai jos: 37 y a A a11 a12 a1n 0 a22 a2n 0 0 a3n 0 0 ann (3.8) Expresia canonică a formei pătratice este: Q(y) 1 y2 1 y2 1 yn (3.9) 11 22 nn 1 a a a unde: y12=a11x1 yn a12x2 a1nxn 22 2 2n n (3.10) a x a x annxn b) Metoda lui Jacobi Fie 1,2,......,n lanţul minorilor diagonali ai matricii A corespunzătoare formei pătratice, daţi de relaţiile: 1 a11 , 2 a11 a21 a12 ,... , n detA (3.11) 22 Dacă toţi i  0 ; i 1,n atunci expresia canonică este dată de formula: Q(y) 1 y2 1 y2 n 1 y2 . (3.12) 1 2 n 1 2 n Observaţii: 1) Expresia canonică a unei forme pătratice nu este unică (coeficienţii i pot avea valori diferite). Totuşi în oricare dintre expresiile canonice asociate aceleiaşi forme pătratice numărul de coeficienţi i pozitivi respectiv negativi este acelaşi; 2) Metoda lui Jacobi nu poate fi aplicată dacă măcar un minor diagonal al matricei A este nul (()i 0); mai mult prin această metodă nu se pot obţine expresiile variabilelor yi definite în (3.7); 38 3) Aplicând metoda lui Jacobi de obţinere a expresiei canonice (3.12), teorema de caracterizare a semnului unei forme pătratice conduce la următoarele concluzii: a) Dacă 1 0,2 0,.........,n 0 atunci Q(x) este pozitiv definită; c) Dacă 1 0,2 0,3 0,......... atunci Q(x) este negativ definită; e) Dacă i  0, ()i 1,n şi în orice altă combinaţie de semne decăt cazul a) sau c) atunci Q(x) este nedefinită ca semn; Evident cazurile b) respectiv d) de semidefinire nu pot fi caracterizate cu metoda lui Jacobi Exemplu. Să se găsească expresia canonică a formei pătratice: Q(x) 2x2 2x1x2 2x2 2x2x3 3x2. 1 2 3 Folosim mai întâi metoda lui Gauss. Obţinem succesiv: 2 1 0 1 1 1 1 1 1 A 1 2 1 1 2 1 0 3 2 0 1 3 0 1 3 0 1 3 : 1 1 1 1 1 1 0 1 4 0 1 4 0 1 3 0 0 7 Expresia canonică este aşadar: Q(y) y2 y2 1 y2 1 2 3 7 Prin metoda lui Jacobi obţinem: 1 2 2 2 1 1 2 2 1 0 3 3 1 2 1 3 0 1 3 Expresia canonică astfel obţinută este deci: Q(z) 1 z2 2z2 3z2 1 2 2 3 3 7 Prin cele două metode am obţinut coeficienţi diferiţi (deci expresii canonice diferite), dar semnele acestora sunt aceleaşi. Forma pătratică este evident pozitiv definită. 39 US4 Programare liniară. Considerații teoretice Principalele noțiuni asociate programării liniare sunt pe larg analizate în această unitate de studiu. După enumerarea (mai mult deducerea) câtorva modele economice remarcabile care conduc la probleme de programare liniară, sunt trecute în revistă formele pe care le poate lua o astfel de problemă. Se acordă atenție și soluțiilor problemelor de programare liniară și proprietăților acestora. O atenție deosebită s-a acordat mulțimilor convexe care constituie o noțiune fundamentală. Principala preocupare a studentului constă în însușirea corectă și completă a noțiunilor teoretice. Puține aplicații concrete sunt accesibile, unele dintre acestea fiind analizate în cadrul activităților tutoriale. Una dintre acestea (programarea în două variabile) care se rezolvă prin metoda grafică este deosebit de importantă deoarece poate servi la exemplificare chestiunilor teoretice introduse. Timpul de lucru este apreciat la șase ore: două ore pentru problema de programare în sine (cu forme diverse de prezentare), două ore pentru prezentarea noțiunii de mulțime convexă și a rezultatelor teoretice asociate și două ore pentru cele trei tipuri de soluții cu proprietățile lor. Orice manual de programare liniară conține chestiunile enumerate aici în diverse maniere. Pozițiile [3] – [9] din lista bibliografică folosesc chiar maniera de prezentare adoptată în lucrarea de față. 40 4.1. MODELE ECONOMICE CE CONDUC LA O PROBLEMĂ DE PROGRAMARE LINIARĂ a) Consumul lunar de materii prime Considerăm un proces de producţie în cursul căruia se realizeză m tipuri de produse P,i 1,m, pentru care se utilizează n tipuri de materii prime Mj, j ,n. Dintr-o unitate de i 1 materii prime Mj se pot produce aij unităţi de produs P . Se doreşte ca numărul de produse de i tipul P (realizate zilnic/săptămânal/lunar /etc.) să fie cel puţin în cantităţile i . Costul unei i b unităţi de materie primă Mj este cj . Se cere să determinăm consumul lunar de materii prime, încât să se realizeze producţia minimă planificată, iar cheltuielile legate de materiile prime să fie minime. În continuare vom elabora modelul matematic al acestui fenomen economic. Notăm cu xj, j 1,n cantităţile de materie primă Mj ce urmează a fi utilizate. Vom scrie sub formă algebrică restricţiile economice impuse mai sus. Aceste restricţii sunt de trei categorii: restricţii legate de planificarea producţiei, restricţii impuse de sensul economic al variabilelor şi restricţia de minimizare a cheltuielilor. Aşa cum aijxj reprezintă cantitatea de produse de tipul i ce se P obţine dintr-o cantitate xj de materie primă Mj , deducem că: i1 1 i2x2 ainxn jn1aijxj ; i 1,n a x a reprezintă cantitatea totală de produse de tipul P , ce se obţine când materiile prime Mj se i utilizează în cantităţile xj, j 1,n . Dar această cantitate nu poate fi mai mică decât i şi deci b obţinem inecuaţiile liniare de forma: j=1aijxjbi bi, i=1,m (4.1) n Cum xj reprezintă mărimi economice (cantităţi de materie primă) rezultă obligatoriu impunerea condiţiilor (de nenegativitate): xj 0, j 1,n (4.2) 41 Cheltuielile totale obţinute ca sumă a cheltuielilor unitare vor fi date de funcţia: f (x) jn1cjxj, x x ,x2, ,xn T n , (4.3) 1 care se cere a fi minimizată. Din punct de vedere matematic problema poate fi formulată astfel: să se găsească în mulţimea soluţiilor sistemului de inecuaţii liniare (4.1) care satisfac condiţiile de nenegativitate (4.2), o soluţie pentru care funcţia f ia valoarea minimă. Problema are sens dacă sistemul (4.1) are cel puţin o soluţie. Inecuaţiile (4.1) le numim restricţiile problemei, (4.2) condiţiile de nenegativitate, iar f poartă denumirea de funcţie obiectiv. b) Folosirea optimă a resurselor Într-un proces de producţie în care se pot realiza n tipuri de articole Aj, j=1,n sunt utilizate m resurse Ri, i=1,m, care sunt limitate de cantităţile bi . Pentru producerea unei unităţi de articol de tipul Aj se consumă o cantitate aij de resursă de tipul Ri . Profiturile unitare pentru articolele Aj presupunem că sunt cj, j=1,n. Să se determine cantitatea xj de articol Aj, j=1,n care urmează a se realiza în procesul de producţie considerat, cu resurse limitate, astfel încât procesul să aibă drept rezultat profit maxim. Limitările resurselor conduc la inecuaţiile: jn1aijxj b , i 1,n (4.4) i Necunoscutele (cantităţile) xj, j 1,n, trebuie să satisfacă condiţiile de nenegativitate: xj 0, j 1,n (4.5) Profitul total este dat de funcţia liniară f (x) jn1cjxj (4.6) care trebuie să ia valoarea maximă. Ne propunem să găsim (măcar) o soluţie x0 x ,x2, ,xn T n (numită în cele ce 1 urmează soluţie optimă) a sistemului de restricţii economice (4.4) care să aibă componentele 42 nenegative, deci să verifice condiţiile (4.5) şi în care funcţia obiectiv (profit) să ia cea mai mare valoare posibilă. Evident, în modelul matematic (corect) al unei probleme economice reale, sistemul de inecuaţii liniare (4.4) are o infinitate de soluţii care verifică condiţiile de nenegativitate (4.5). Dorim să găsim acea soluţie x0 pentru care are loc: (min) f (x) f ( x0 ) b) Problema amestecului optim Să presupunem că pentru realizarea unui amestec se pot întrebuinţa n materii prime Mj, j 1,n, disponibile în candităţi nelimitate, ce conţin m substanţe S , i 1,m. O unitate de i materie primă Mj , conţine cantitatea aij din substanţă Si . Amestecul trebuie să conţină substanţa Si în exact cantitatea b . Costurile unitare ale materiilor prime Mj sunt ci unităţi i băneşti. Să de determine cantităţile xj, j 1,n de materie primă Mj necesare realizării amestecului cu compoziţia prescrisă şi la un cost total minim. Restricţiile sunt şi acum de trei feluri: restricţiile de compoziţie a amestecului, restricţiile datorate sensului economic al variabilelor şi restricţia de minimizare a costului total. Traducerea algebrică a primelor două categorii este: jn1 ijxj b , i 1,m (4.7) a i xj 0, j 1,n (4.8) Costul total al materiilor prime necesare pentru realizarea unei unităţi de amestec este dat de funcţia: f (x) jn1cjxj (4.9) care se cere a fi minimizată. 43 4.2. FORMULAREA GENERALĂ A UNEI PROBLEME DE PROGRAMARE LINIARĂ Se constată că modelele matematice particulare, ataşate problemelor economice prezentate anterior prezintă o serie de asemănări, fapt ce permite înglobarea lor într-un model general. Să considerăm un sistem care poate conţine ecuaţii, inecuaţii de un semn sau altul, adică un sistem de forma: 11 1 12x2 1nxn 1 21 1 22 2 2n n a x a a b a x a x a x b .............................................. 2 (4.10) am1 x am2x2 amnxn bm 1 numit sistemul de restricţii (economice) al problemei, care poate fi scris condensat (utilizând indicele de sumare) şi sub forma: jn1 ijxj b , i 1,m (4.11) a i (semnul înseamnă semnul =, ori , ori , ori <, ori >). Condiţiile de tipul: xj 0, j 1,n (4.12) sunt numite condiţii de nenegativitate şi sunt necesare datorită caracterului economic real al variabilelor (d.p.d.v. strict matematic putem rezolva şi problemele în care aceste condiţii nu apar). Funcţia liniară: f (x) jn1cjxj (4.13) care trebuie să o optimizăm (minimizăm sau maximizăm) se numeşte funcţie obiectiv (funcţie cost sau funcţie profit conform tipului de problemă studiat). Problema generală a programării liniare se formulează astfel: în mulţimea soluţiilor sistemului (4.11) ce satisfac condiţiile (4.12) să se găsească o soluţie (numită soluţie optimă) ce dă valoarea minimă sau maximă funcţiei definite de (4.13). 44 1 a P a n 1 , a 4.3. FORME ALE UNEI PROBLEME DE PROGRAMARE LINIARĂ Modelul matematic al unei probleme de programare liniară poate fi scris sub diverse forme potrivit scopurilor particulare urmărite. Aceste forme au doar o semnificaţie strict matematică şi uşurează respectiv simplifică demonstraţiile anumitor rezultate matematice (propoziţii sau teoreme) absolut necesare în fundamentarea modului (algoritmului) de rezolvare a problemelor de programare liniară. a) Forma matriceală Considerând matricele: A (aj )i1,m Mm,n , x (x ,x2,.....,xn)T  n , j1,n i 1 b(b ,b ,....,bm)T  m şi C (c ,c2,....,cn) n , putem scrie problema de programare liniară 1 2 1 sub forma: (min/max) f(x) = C x Ax b , (4.14) x 0 b) Forma vectorială Dacă introducem vectorii: P 11,a21, ,am1 T , 2 12,a22, ,am2 T , ,P a n,a2n, T mn P0 = b. atunci problema poate fi scrisă sub forma: (min/max) f(x) = C x x P0 x2P xnP P (4.15) x 1 1 2 n 0 c) forma canonică O problemă de programare liniară este sub formă canonică dacă ea se scrie: (min) f(x) = Cx Ax b (4.16) x 0 45 sau (max) f(x) = Cx Ax b (4.17) x 0 d) forma standard O problemă de programare liniară este sub formă standard dacă restricţiile ei sunt sub formă de egalităţi. O problemă de programare liniară poate fi scrisă sub forma standard matricială sau vectorială, adică: (min/max) f(x)=cx Ax b sau x 0 (min/max) f(x) = C x x P0 x2P xnP P (4.18) x 1 1 2 n 0 Observaţie. Formele (4.16), (4.17), (4.18) ale unei probleme de programare liniară sunt echivalente dacă ţinem seama de următoarele: 1. Sensul unei inegalităţi se schimbă prin înmulţirea cu –1; 2. O inecuaţie de forma i1 1 i2x2 inxn i , poate fi scrisă ca o ecuaţie a x a a b i1 1 i2x2 inxn xc b , prin adăugarea a unei noi variabile (variabilă de a x a a i i compensare) xc 0; i 3. O inecuaţie de forma: i1 1 i2x2 inxn i , poate fi scrisă ca o ecuaţie a x a a b i1 1 i2x2 inxn xc b , scăzând o variabilă de compensare xc 0 ; a x a a i i i 4. O ecuaţie poate fi scrisă ca două inecuaţii de semn contrar ; 5. Problemele de maxim pentru funcţia f (x) pot fi reduse la o problemă de minim pentru funcţia f '(x) f (x); cu alte cuvinte în cazul unei funcţii (forme) liniare de tipul (4.13) are loc egalitatea : (max)f (x)  (min)(f (x))  (min)f '(x) , pentru ()x n Având în vedere cele de mai sus, în continuare ne vom ocupa numai de probleme aflate sub forma standard cu cerinţa de minim pentru funcţia obiectiv. 46 i i i i 4.4. MULŢIMI CONVEXE Fie spaţiul vectorial n şi M o submulţime nevidă a sa (M n , M ). Definiţia 4.1. Mulţimea M n se numeşte convexă dacă ( )x1 ,x2 M şi ( )λ 0,1 avem: x λx1 (1 λ)x2 M (4.19) Definiţia 4.2. Un punct (vector) x M se numeşte punct extrem sau vârf al mulţimii convexe M, dacă nu există x1 x2 din M şi nu există λ (0,1) astfel încât x λx1 (1 λ)x2 (4.20) în caz contrar, punctul se numeşte punct interior. Definiţia 4.3. Numim combinaţie liniară convexă a vectorilor xi M, i=1,k o expresie de forma: k λ x , cu λ i 1 [0,1] k λ 1 (4.21) i=1 Observaţie. Se poate demonstra (relativ uşor) că orice punct interior se poate scrie (exprima) întotdeauna ca o combinaţie liniară convexă de punctele extreme ale mulţimii convexe; evident, conform definiţiei 4.2, un vârf (punct extrem) al unei mulţimi convexe nu poate fi scris ca o combinaţie liniar convexă de puncte interioare mulţimii. 4.5. SOLUŢII ALE UNEI PROBLEME DE PROGRAMARE LINIARĂ Fie o problemă de programare liniară aflată sub forma standard matricială: (min) f(x)=cx Ax b (4.22) x 0 47 1 2 Definiţia 4.4. Un vector x n ale cărui componente xj, j ,n verifică sistemul de 1 restricţii şi condiţiile de nenegativitate se numeşte soluţia admisibilă a unei probleme de programare liniară. Mulţimea soluţiilor admisibile o vom nota cu Sa . Teorema 4.1. Mulţimea soluţiilor admisibile Sa este o mulţime convexă. Demonstraţie. Fie x1 x2 Sa două soluţii admisibile. Atunci avem: Ax1 b, x1 0 ; Ax2 b, x2 0. Să considerăm un vector x de forma: x λx1 (1 λ)x2, cu λ 0,1 . Să arătăm că x Sa . Într-adevăr vom putea scrie: Ax A λx1 (1 λ)x2 = λAx1+(1 λ)Ax2= λb+(1 λ)b = b şi deci verifică sistemul de restricţii. Cum λ 0,1 , iar x1, x2 0, rezultă λx1 (1 λ)x2 = x 0 şi deci x este soluţie admisibilă a problemei de programare liniară. Definiţia 4.5. O soluţie admisibilă care are cel mult m componente strict pozitive, iar vectorii formaţi cu coloanele matricei A, corespunzătoare acestor componente sunt liniar independenţi, se numeşte soluţie admisibilă de bază. Dacă soluţia admisibilă de bază are exact m componente strict pozitive, spunem că este nedegenerată, iar dacă are mai puţine, degenerată. Mulţimea soluţiilor admisibile de bază o vom nota cu Sab şi este o mulţime finită (card Sab Cn ) deci evident, nu este o mulţime convexă. m Exemplul 4.1. Fie problema de programare: (min) f(x) 2x x 2x xx2 2x3 xj 0, j 1,4 1 1 2 3 2 x x 2x4 2 x4 4 48 În baza definiţiilor (4.4) şi (4.5) constatăm că: x1= 2/3, 1/3, 2/3, 5/3 T este soluţie admisibilă x2= 0, 1, 0, 5 T este soluţie admisibilă de bază nedegenerată x3= 0, 0, 2, 0 T este soluţie admisibilă de bază degenerată x4= 1, 0, 1, 0 T este soluţie admisibilă, dar nu este de bază deoarece (vectorii P1 şi P3 nu sunt liniar independenţi) Teorema 4.2. Orice soluţie admisibilă de bază a unei probleme de programare liniară, este un punct extrem al mulţimii soluţiilor admisibile Sa şi reciproc, orice punct extrem al mulţimii soluţiilor admisibile este o soluţie de bază a problemei de programare liniară, cu alte cuvinte are loc: xSab  xSa este punct extrem (4.23) Definiţia 4.6. O soluţie admisibilă x se numeşte soluţie optimă, dacă realizează minimul sau maximul funcţiei obiectiv f; mulţimea soluţiilor optime o vom nota cu SO Se poate demonstra că mulţimea soluţiilor optime este o mulţime convexă. Evident, teoretic SO poate să fie vidă (problema de programare aferentă nu are soluţii optime), să aibă o infinitate de elemente (în acest caz spunem că problema de programare aferentă are o infinitate de soluţii optime). Dar cazul cel mai des întâlnit şi cel care apare în imensa majoritate a problemelor economice reale este cel în care SO are numai un element (problema de programare aferentă are soluţie optimă unică). Teorema 4.3. Dacă o problemă de programare are o soluţie optimă finită, atunci există un punct extrem al mulţimii soluţiilor admisibile Sa în care funcţia obiectiv ia aceeaşi valoare, adică: (min)f (x)  fO  (finit)  ()xSab a.i. f (x)  fO (4.24) 49 US5 Metoda Simplex În unitatea de studiu precedentă s-a prezentat modelul de bază al unei probleme de programare liniară. În unitatea prezentă dăm o metodă de rezolvare a acesteia: metoda Simplex. Aici se vede importanța unităților de studiu anterioare (US1 și US2) deoarece se lucrează cu vectori (elemente ale spațiului vectorial al restricțiilor) și folosim (pentru determinarea componentelor soluției) metoda transformărilor elementare. De mare utilitate este algoritmul prezentat în US2 care se aplică aici fără modificare. Se dau principiile algoritmului Simplex (criteriile de optim, de intrare, de iesire etc.) și se descriu etapele algoritmului. Se prezintă apoi cazul în care soluția inițială trebuie construită (metoda celor doua faze). Studenții trebuie să rețină mai întâi criteriile, apoi au de reținut forma tabelului pe care se aplică algoritmul Simplex, cu precizarea clară a modului de completare a fiecărui element al acestuia. Ulterior ei trebuie să deprindă ordinea în care se aplică cele trei criterii fundamentale și calculele care se fac pentru fiecare. Ultimul lucru îl constituie efectuarea calculeleor (transformările elementare) pentru obținerea unei soluții mai bune. După fixarea temeinică a algoritmului se poate trece la metoda celor două faze. Sunt necesare zece ore pentru parcurgerea acestui material (două ore pentru parcurgerea principiilor, o oră pentru descrierea algoritmului, patru ore pentru aplicații ale algoritmului, trei ore pentru metoda celor două faze). Din lista bibliografică sunt utile pozițiile [3] – [6], [8], [9]. În [7] este prezentată o variantă puțin modificată dar interesantă. 50 5.1. PRINCIPIILE METODEI SIMPLEX Să considerăm problema de programare liniară sub forma standard matricială: (min) f(x)=cx (*) Ax b (**) (5.1) x 0 (***) sau vectorială (min) f(x) = C x x P 0 x2P xnP 1 1 2 x n (*) P (**) (5.2) (***) 0 Raţionând în mod simplist, pentru a rezolva o problemă de programare liniară (adusă la forma standard), ar trebui să parcurgem următorii paşi: 1. Să determinăm toate soluţiile sistemului liniar compatibil nedeterminat (**), lucru relativ uşor şi cunoscut încă din liceu; 2. Să eliminăm (cum!?) din infinitatea de soluţii a lui (**) acele soluţii (tot o infinitate) care nu verifică condiţiile de nenegativitate (***); 3. Din infinitatea de soluţii admisibile rămase (adică soluţiile sistemului (**) care satisfac condiţiile (***)) să o găsim (cum!?) pe aceea sau acelea care fac funcţia obiectiv să ia cea mai mică valoare, adică care satisfac(e) şi condiţia (*). Rezultatele indicate în capitolul precedent ne oferă un mod simplu şi elegant de determinare a soluţiei optime a unei probleme de programare liniară. Astfel, teoremele 3.1.2 şi 3.1.3 ne indică că trebuie să căutăm soluţiile optime ale unei probleme de programare liniară, în mulţimea soluţiilor admisibile de bază (Sab este o mulţime finită având cel mult Cm elemente). Cu alte cuvinte, va trebui să parcurgem următorii paşi: n 1. Determinăm (eventual folosind transformări elementare) toate soluţiile de bază ale sistemului (**); sunt cel mult Cm soluţii de bază, deci un număr finit (şi nu infinit); n 2. Eliminând soluţiile cu componente negative (neadmisibile), care nu satisfac condiţiile (***), obţinem mulţimeaSab a soluţiilor de bază admisibile, printre care se află şi soluţi- a/ile optim-ă/e căutat-ă/e; 51 3. Calculăm valoarea funcţiei obiectiv în toate soluţiile admisibile de bază aflate la pct. 2). Soluţia optimă va fi cea în care funcţia obiectiv f (x)ia valoarea minimă. Totuşi, chiar şi în cazul problemelor în care numărul restricţiilor (ecuaţiilor)m respectiv al variabilelor (necunoscutelor) n este relativ mic, numărul de calcule necesar aplicării metodei de mai sus este extrem de mare, ceea ce o face aproape imposibil de aplicat în realitate. De exemplu, pentru m=10, n=40 numărul de soluţii de bază al sistemului (**), care trebuie determinate, este de C10 846.660.528 soluţii (!!). 40 Vom prezenta în continuare o metodă (algoritm) numită metoda (algoritmul) SIMPLEX, elaborată de matematicianul german G.B. Dantzig în anul 1906, algoritm care are avantajul unui volum extrem de redus de calcule şi utilizarea numai a transformărilor elementare ca aparat matematic. Să presupunem fără a restrânge generalitatea că matricea lărgită a sistemului de restricţii este de forma: 1 0 0 1m+1 1n 1 0 1 0 2m+1 2n 2 (5.3) 0 0 1 mm+1 mn m Vectorii P1, P2, . . ., Pm formaţi cu primele m coloane ale matricei A, fiind liniar independenţi, formează o bază în m . În acest caz un vector oarecare P , j=m+1,n se exprimă j în mod unic cu ajutorul vectorilor bazei sub forma: Pj=α1jP +α2jP + +αmjPm i 1αij i (5.4) 1 2 m P Să notăm cu CB vectorul linie format cu primele m componente ale vectorului C, adică coeficienţii variabilelor x ,x2, ,xn din funcţia obiectiv ce corespund vectorilor care au intrat în 1 bază. Introducem mărimile: zj-cj=zj-CBP i=1ciαij ; j=1,n (5.5) j m Dacă presupunem acum că βi, i=1,m sunt pozitive şi dacă avem în vedere (5.3) deducem că vectorul x de componente, 52   x1 = β1,x2 = β2, ,xm= βm, xm+1= 0,…, xn = 0, este o soluţie admisibilă de bază a problemei de programare considerată. Teorema 5.1. (Criteriul de optim) Dacă diferenţele zj cj 0, j = m+1,n, atunci problema de programare liniară are optim finit şi x este soluţia optimă căutată. Observaţii: a) se constată uşor că zj cj  0 pentru j1,m b) În problemele de programare, în care avem cerinţa de maxim pentru funcţia obiectiv, criteriul de optimalitate va fi: zj cj  0, j m1,n c) Din demonstraţia teoremei rezultă că în cazul în care există zj cj  0 în probleme de minim pentru funcţia obiectiv pot exista soluţii admisibile pentru care funcţia obiectiv să ia o valoare mai mică decât pentru x, deci x nu este optimă. În cazul în care x nu este soluţie optimă vom căta să obţinem din vechea soluţie de bază admisibilă x o nouă soluţie de bază admisibilă y “mai bună” decât vechea soluţie x, în sensul că valoarea funcţiei obiectiv f (x)este mai “bună” (adică mai mică) decât în vechea soluţie ( f (y)  f (x)). Dantzig demonstrează că, acest lucru se poate obţine înlocuind un vector Pi (corespunzător unei variabile principale a soluţiei x)din bază, cu un vector Pj (corespunzător unei variabile secundare a soluţiei x) din afara bazei. Să notăm cu yk , k 1,2, ,i 1, j,i 1, ,m componenetele noii soluţii de bază y . Din (1.2.5) înlocuind bk cu k , aki cu ij , xjk cu yk , xj cu j obţinem: yk  k  i kj , ij  0, k 1,m ; k  i ij yj  i , k i ij      (5.6) 53 Vectorul Pi (care iese din bază) şi vectorul Pj (care va intra în bază) se determină pe baza următoarelor două teoreme: Teorema 5.2. (Criteriul de ieşire din bază). Condiţiile necesare pentru înlocuirea vectorului Pi cu vectorul Pj sunt: i  i  mink  k , k 1,m ; kj  0, (5.7) ij kj           Să vedem acum ce condiţii trebuie să îndeplinească un vector Pj pentru a intra în bază în locul vectorului Pi pentru a obţine o soluţie "mai bună". Teorema 5.3. (Criteriul de intrare în bază) Condiţia pe care trebuie să o îndeplinească un vector Pj, j m1,n pentru a intra în bază în locul vectorului Pi în vederea obţinerii unei soluţii "mai bune" este: zj cj  0 (5.8) 5.2. ETAPELE ALGORITMULUI SIMPLEX PRIMAL Rezultatele enunţate în capitolul anterior se concretizează în următorul mod de lucru, cunoscut în literatura de specialitate sub denumirea de algoritmul Simplex primal. Acesta constă dintr-o succesiune de iteraţii, care îmbunătăţesc treptat o soluţie admisibilă de bază iniţială. Pe baza consideraţiilor din aliniatul 5.1. putem pune în evidenţă etapele algoritmului care ne vor conduce la soluţia optimă a problemei dacă aceasta există. 1. Se determină o soluţie admisibilă de bază, fie aceasta x(oare e atât de simplu!?); 2. Se construieşte tabelul Simplex corespunzător acestei soluţii (este redat mai jos); 3. Se calculează diferenţele zj cj , j1,n aflate pe ultima linie a tabelului iniţial. Aplicăm criteriul de optim. Putem avea următoarele patru situaţii: a) dacă zj cj  0, j unică!); b) dacă zj cj  0, j m 1,n atunci soluţia admisibilă de bază x este optimă (şi m 1,n şi ()zj cj  0iar vectorul j corespunzător are şi P componente strict pozitive, atunci soluţia admisibilă de bază x este optimă (dar nu este unică; problema admite o infinitate de soluţii optime cu aceeaşi valoare a funcţiei obiectiv); 54 k  c) Dacă există indici j, j m 1,n pentru care zj cj  0, iar toate componentele vectorului Pj corespunzător sunt mai mici sau egale cu zero (ij  0), atunci problema de programare are optim infinit; d) Dacă ()j ; j m 1,n pentru care zj cj  0 şi nu toţi ij  0, atunci se trece la etapa următoare (soluţia nu este optimă). 4. Se aplică criteriul de intrare, determinând vectorul Pj care intră în bază ca fiind acel vector pentru care: zj cj max (z zk ck 0 ck) 5. Se aplică criteriul de ieşire. Părăseşte baza vectorul Pi pentru care i min k ij kj 0 kj  i  min k 0 6. Se face schimbarea de bază obţinând o nouă soluţie admisibilă de bază y , pentru care etapele algoritmului se reiau. Datele problemei de programare liniară şi elementele ce au intervenit în teoremele 5.1, 5.2, 5.3 respectiv în algoritmul prezentat mai sus, le sintetizăm în următorul tabel Simplex: c1 B CB P0 P 1 P1 c1 1 1 P2 c2 2 0 ⋮ ⋮ ⋮ ⋮ Pi ci i 0 ⋮ ⋮ ⋮ ⋮ Pm cm m 0 zj-cj - 0 c2 ... ci P2 ..... Pi 0 ⋯ 0 1 ⋯ 0 ⋮ ⋮ 0 ⋯ 1 ⋮ 0 ⋯ 0 0 ... 0 ... ... . ⋯ ⋯ ⋯ ⋯ ⋯ cm Pm 0 0 ⋮ 0 ⋮ 1 0 ... ... . ⋯ ⋯ ⋯ ⋯ ... cj ... cn Pj ... Pn α1j ... α1n α2j ... α2n ⋮ ⋮ αij ... αin ⋮ ⋮ αmj ... αmn zj-cj ... zn-cn k  O j P P 1 2 i m 55 j i ij Observaţii: a) Valoarea funcţiei obiectiv pentru soluţia admisibilă de bază este dată de: f(x) = i=1ciβi = CBP m 0 b) Diferenţele zj cj se calculează cu relaţiile: zj c = m c α i=1 cj i) Dacă soluţia nu este optimă (există zj cj  0) alegem diferenţa pozitivă care este cea mai mare, determinând în felul acesta vectorul Pj care intră în bază. ii) Dacă există mai multe astfel de diferenţe, atunci poate intra în bază oricare dintre vectorii ce corespund acestora. iii) Dacă pe coloana vectorului Pj care trebuie să intre în bază toate elementele sunt mai mici sau egale cu zero, atunci problema admite optim infinit. iv) Dacă există şi elemente strict mai mari ca zero considerăm toate rapoartele dintre componentele lui P0 şi componentele strict pozitive ale vectorului Pj, rapoarte de forma i i ij componentele lui P componentele lui P O j 0 cu αij>0. Cel mai mic raport i  0 va indica vectorul care părăseşte baza. Dacă vectorul Pi părăseşte baza atunci noua soluţie admisibilă de bază se obţine cu ajutorul unui pivotaj cu elementul pivot αij. Exemplu 5.1. Pentru fabricarea a 3 tipuri de produse Pj, j=1,3 se utilizeată 2 tipuri de resurse R1 şi R2. Consumurile specifice, cantităţile de resurse şi profiturile unitare sunt date în tabelul: P R P1 P2 R1 1 1 R2 2 3 Profit 4 1 Limitări P3 resurse 2 30 1 57 5 56 1 2 1 2 3 1 1 c Să se determine cantităţile de produse Pj, j=1,3 ce urmează a fi realizate pentru a obţine un profit maxim. Soluţie. Notăm cu x ,x2,x3 cantităţile de produse de tipul P1, P2, P3 ce urmează a fi 1 realizate. Modelul matematic al problemei este: (max) f(x) 4x x x x 2x 30 2x 3x2 x3 57 xj 0, j=1,3 5x3 Aducem problema la forma standard folosind variabile de compensare xc şi xc şi 4 5 considerăm problema de minim pentru funcţia –f . Vom avea problema: (min) ( f(x)) 4x x2 5x3 x x2 2x3 x4 30 2x 3x2 x3 x5 57 1 c 1 c c xj 0, j 1,3, x4 0, x5 0 O soluţie admisibilă de bază este x = 0, 0, 0, 30,57 T . Tabelul simplex se prezintă astfel: B CB P0 Pc 0 30 4 Pc 0 57 5 zj-cj – 0 P3 – 5 15 Pc 0 42 zj-cj – –75 5 P3 – 5 1 P1 – 4 28 zj-cj – –117 – 4 – 1 – 5 P1 P2 P3 1 1 ② 2 3 1 4 1 5 1/2 1/2 1 3/2 5/2 0 3/2 – 3/2 0 0 – 1/3 1 1 5/3 0 0 –1 0 0 0 c c 4 5 P P 1 0 0 1 0 0 1/2 0 – 1/2 1 – 5/2 0 2/3 –1/3 –1/3 2/3 –2 –1 i 1  30 15 2 2  57 57 1 30 28 - - - 57 1 1 2 3 4 1 2 În prima etapă au intrat în bază vectorii xc şi xc . Pe linia diferenţelor zj cj observăm 4 5 că există diferenţe strict pozitive, deci conform cazului 3)d de la criteriul de optim soluţia iniţială nu este optimă. Cea mai mare diferenţă pozitivă este z3-c3 = 5, şi conform criteriului de intrare, vectorul P3 va intra în bază. Avem rapoartele 4 4 43 30 2 15, 5 5 53 57 1 57. Cel mai mic raport este 4 15, ceea ce înseamnă (conform criteriului de ieşire)că vectorul Pc părăseşte baza fiind 4 înlocuit de vectorul P3 . Un pivotaj cu element pivot 2 (încercuit) conduce la o nouă soluţie admisibilă de bază x 0, 0, 15, 0, 42 T care nu este optimă (există diferenţe zj cj  0 în tabelul corespunzător noii soluţii găsite). Continuând cu o nouă iteraţie cu elementul pivot 3/2 (încercuit), obţinem soluţia admisibilă de bază x 28, 0, 1, 0, 0 T care este optimă şi unică. Lăsând la o parte variabilele de compensare introduse, obţinem soluţia optimă a problemei iniţiale: x 28, 0, 1 T . Din punct de vedere economic aceasta înseamnă realizarea a 28 unităţi de produs P1 şi a unei unităţi de produs P3 , profitul maxim posibil de realizat fiind de 117 (u.m.). Exemplul 5.2 Să se determine soluţiile optime ale problemei: (min) f(x) 3x 2x x 2x x x 7, x 4x2 x3 x5 1 xj 0, j 1,5 4x3 O soluţie admisibilă de bază iniţială este x 0, 0, 0, 7, 1 T . Tabelul Simplex corespunzător este: 58 B CB P0 P 0 7 4 P 0 1 5 zj-cj – 0 P 0 8 4 P – 4 1 3 zj-cj – – 4 – 3 – 2 – 4 P1 P2 P3 1 2 – 1 1 – 4 ① 3 2 4 2 – 2 0 1 – 4 1 –1 18 0 0 0 c c 4 5 P P 1 0 0 1 0 0 1 1 0 1 0 – 4 i 4  0 5 1 4  0 3  0 Algoritmul se opreşte, chiar dacă soluţia admisibilă de bază găsită x 0, 0, 1, 8, 0 T nu este optimă (diferenţa z2 c2 180), deoarece vectorul P2 care ar trebui să intre în bază are toate componentele negative. Problema de programare liniară considerată admite optim infinit conform cazului 3)c. 5.3. METODA CELOR DOUĂ FAZE Algoritmul Simplex este un algoritm rapid convergent la soluţia optimă. Astfel, plecând de la o soluţie de bază admisibilă iniţială a sistemului în formă standard, cu m ecuaţii şi n necunoscute, după cel mult m+n iteraţii obţinem soluţia optimă. De exemplu, în cazul unui sistem cu 10 ecuaţii şi 40 de necunoscute, plecând de la o soluţie admisibilă iniţială, în cel mult 10+40=50 de paşi (iteraţii) se ajunge la soluţia optimă (dacă există). Determinarea (eventual prin transformări elementare) a unei soluţii de bază pentru un sistem liniar, chiar de dimensiuni mari, este extrem de simplă. Problema care apare, este admisibilitatea acesteia. De exemplu, în cazul sistemului prezentat anterior cu C10 846.660.528 40 soluţii de bază, se poate întâmpla ca foarte multe dintre acestea, sau chiar toate, să fie neadmisibile. Am putea astfel calcula un număr imens ( sute de milioane) de soluţii neadmisibile (deci inutile), până obţinem acea soluţie admisibilă iniţială. După găsirea acesteia, aplicând algoritmul Simplex vom obţine soluţia optimă foarte rapid, în cel mult 50 de iteraţii. 59 Vom indica o metodă mai lesne de aplicat în practică pentru a obţine o soluţie admisibilă de bază iniţială: metoda celor două faze. Să considerăm că problema de programare liniară (5.1) este sub formă standard, termenii liberi sunt pozitivi şi că matricea A a sistemului de restricţii nu conţine nici o coloană a matricei unitate. Vom adăuga la fiecare restricţie o variabilă nenegativă, numită variabilă artificială. Fiind m restricţii vom adăuga m variabile artificiale xm , k=n+1, n+m , şi considerăm o nouă problemă k de programare liniară, numită problema artificială ataşată problemei iniţiale, de forma: (min) g(xa )= n+m xa k k=n+1 a Ax + Imx = b (5.9) x 0,x 0 a Aici Im este matricea unitate de ordinul m şi xa xn 1,xn 2,.......xn m T vectorul variabilelor a a a artificiale. Se constată că orice soluţie admisibilă a problemei (5.9) în care variabilele artificiale sunt nule, devine după înlăturarea acestora o soluţie admisibilă de bază a problemei (5.1). Metoda celor două faze constă în: Faza întâi. Se determină soluţia optimă a problemei (5.9) folosind algoritmul Simplex. Soluţia de bază admisibilă iniţială artificială, se determină direct din sistemul artificial, considerând variabile principale pe cele artificiale. Aplicând problemei (5.19) algoritmul Simplex primal, ajungem la una din situaţiile: a) (min) g(xa )=0 şi nici un vector corespunzător variabilelor artificiale nu este în bază. b) (min) g(xa )=0 şi în bază rămân şi vectori corespunzători variabilelor artificiale. c) (min) g(xa ) 0. d) Faza a doua. În situaţia a) soluţia optimă a problemei (5.9), după înlăturarea variabilelor artificiale, reprezintă o soluţie de bază admisibilă, iniţială, nedegenerată pentru problema (5.1). Cu această soluţie iniţială, se aplică algoritmul Simplex primal în vederea găsirii soluţiei optime a problemei iniţiale. Din tabelul final al fazei întâi se elimină elementele ce ţin de variabilele 60 artificiale (coloanele vectorilor artficiali), se modifică coeficienţii funcţiei obiectiv şi se calculează noile diferenţe zj cj , obţinându-se primul tabel Simplex al fazei a doua. . Situaţia b) conduce la o soluţie admisibilă de bază degenerată pentru problema iniţială. Consideraţiile anterioare rămân valabile, însă nu vor fi eliminate coloanele corespunzătoare variabilelor artificiale, pentru care vectorii artificiali corespunzători au rămas în bază. În situaţia c) problema iniţială nu admite soluţii admisibile. Exemplul 5.3. Să se determine soluţiile optime ale problemei de programare liniară: x 3x 2x 3, 3x 3x2 x3 4 1 2 3 1 x 0, j=1,3 j (min) f(x) 3x 4x2 2x3 1 Aducem problema la forma standard, ea devenind: x 3x2 2x3 x4 3, 3x 3x2 x3 x5 4, c 1 c 1 c c j 4 5 x 0, j=1,3, x 0, x 0 (min) f(x) 3x 4x2 2x3 1 Vom întrebuinţa metoda celor două faze. În prima fază căutăm soluţia optimă a problemei (artificiale): x 3x2 2x3 x4 x6 3, 3x 3x2 x3 x5 x7 4, c a 1 c a 1 c c a a xj 0, x4, x5, x6 , x7 0 (min) g(xa ) xa xa 6 7 Soluţia admisibilă de bază de pornire este x 0, 0, 0, 0, 0, 3, 4 T . 61 Tabelul simplex are forma: B CB P0 Pa 1 3 6 Pa 1 4 7 zj-cj – 7 P2 0 1 Pa 1 1 7 zj-cj – 1 P2 0 5/6 P1 0 1/2 zj-cj – 0 0 0 P1 P2 1 ③ 3 3 4 6 1/3 1 ② 0 2 0 0 1 1 0 0 0 0 P3 2 1 3 2/3 –1 –1 5/6 – 1/2 0 0 0 c c 4 5 P P –1 0 0 –1 –1 –1 –1/3 0 1 –1 1 –1 –1/2 1/6 1/2 –1/2 0 0 1 1 a a 6 7 P P 1 0 0 1 0 0 1/3 0 –1 1 –1 0 1/2 –1/6 –1/2 1/2 –1 –1 Am găsit soluţia optimă x 1/2, 5/6, 0, 0, 0, 0, 0 T a problemei artificiale. Soluţia admisibilă de bază pentru problema iniţială va fi: x 1/2, 5/6, 0, 0, 0 T , iar sistemul de restricţii cel ce rezultă din ultima etapă a algoritmului simplex pentru problema artificială. Tabelul Simplex corespunzător fazei a doua are forma: B CB P2 –4 P1 3 zj-cj – P2 –4 Pc 0 4 zj-cj – P0 5/6 1/2 – 11/6 4/3 1 16/3 3 – 4 2 P1 P2 P3 0 1 5/6 1 0 –1/2 0 0 –41/6 1 1 1/3 2 0 –1 –7 0 –10/3 0 c 4 P –1/2 1/2 7/2 0 1 0 0 c 5 P 1/6 –1/2 – 13/6 –1/3 –1 4/3 62 Problema admite optim infinit deoarece vectorul Pc nu are şi coordonate strict pozitive. 5 Exemplul 5.4. Să se determine soluţiile optime ale problemei de programare liniară: x 2x x 2x 4 3x x x 8x 3 1 4 5 6 2 4 5 6 x2 x3 x4 x5 3x6 2 xj 0, j 1,6 (min)f(x) 2x1 3x2 x3 4x4 5x5 8x6 Prima şi a treia coloană a matricei A sunt coloane ale matricei unitate de ordinul trei şi atunci vom adăuga la restricţia a doua o variabilă artificială xa . 7 Vom căuta în prima fază soluţia optimă a problemei: x1 2x4 x5 2x6 4, a 2 4 5 6 7 3x x x 8x x 3, x2 x3 x4 x5 3x6 2. a xj 0, j 1,6 , x7 0 (min)g(xa ) xa 7 O soluţie admisibilă de bază pentru această problemă este: x 4, 0, 2, 0, 0, 0, 3 T Tabelul simplex are forma: B CB P0 P1 0 4 Pa 1 3 7 P3 0 2 zj-cj – 3 P1 0 4 P2 0 1 P3 0 1 zj-cj – 0 0 0 0 0 P1 P2 P3 P4 1 0 0 –2 0  0 –1 0 1 1 1 0 3 0 –1 1 0 0 –2 0 1 0 –1/3 0 0 1 4/3 0 0 0 0 0 0 P5 P6 1 2 1 –8 –1 –3 1 –8 1 2 1/3 –8/3 –4/3 –1/3 0 0 1 a 7 P 0 1 0 0 1 1/3 –1/3 -1 63 Soluţia admisibilă de bază pentru problema iniţială va fi x 4, 1, 1, 0, 0, 0 T . Tabelul simplex pentru problema iniţială (faza a doua) are forma: B CB P0 P1 –2 4 P2 –3 1 P3 –1 1 zj-cj – –12 –2 –3 –1 4 P1 P2 P3 P4 1 0 0 –2 0 1 0 –1/3 0 0 1 4/3 0 0 0 –1/3 5 8 P5 P6 1 2 1/3 –8/3 –4/3 –1/3 –20/3 –11/3 Soluţia optimă a problemei este xopt (4,1,1,0,0,0)T , (min)f(x) = –12. 64 US6 Probleme de transport Problemele de transport sunt probleme de programare liniară particulare. Ele se referă la distribuirea unui produs P aflat în cantități date într-un număr de depozite și solicitat de un număr de centre de desfacere, de asemenea în cantități cunoscute. Această problemă nu se rezolvă cu algoritmul Simplex din cauza dimensiunilor matricii asociate. Algoritmul prezentat are aceleași caracteristici cu algoritmulk Simplex: un criteriu de optim, un criteriu de intrare și un criteriu de ieșire. Toate acestea sunt descrise în prezenta unitate de studiu. Este descris apoi algoritmul, modul în care se obține o soluție mai bună dintr-o soluție dată. Studentul trebuie să învețe mai întâi să echilibreze problema. Ulterior el trebuie sa obțină o soluție nedegenerată printr-una dintre cele două metode uzuale (metoda diagonalei și metoda costului minim). După verificarea optimalității soluției, urmează construirea unei soluții mai bune. Calculele nu sunt complicate, dar algoritmul pune în general probleme studenților. Sunt necesare două ore pentru însușirea noțiunilor teoretice și încă patru ore pentru aplicarea cursivă și ușoară a criteriilor. Titlurile bibliografice [3] – [9] conțin tratarea subiectului în manera pe care și noi am abordat-o. * * * Să considerăm că un anumit produs este stocat în m depozite Di în cantităţile ai , i 1,m. Produsul este solicitat de n centre de consum Ci în cantităţile bj , j1,n . Costul transportului unei unităţi de produs de la depozitul Di la centrul Cj este aij u.m. Se cere a se determina cantităţile xij de produs care urmează a fi transportate de la depozite la centrele de consum astfel ca disponibilul a fie epuizat, în fiecare depozit, cererea să fie satisfăcută exact, în fiecare centru de consum, iar costul total al transportului produsului să fie minim. Se presupune că m a  n b i1 1   i j , adică disponibilul din depozite este egal cu cererea totală a centrelor de consum. O astfel de problemă este denumită problema de transport echilibrată. 65 x , Restricţiile problemei sunt din nou de trei tipuri: restricţii asupra disponibilului şi cererii, restricţiile datorate sensului economic al variabilelor şi restricţia minimizare a costului. Datele modelului economic le prezentăm într-un tabel de forma: C D C1 C2 … Cj … Cn ai D1 D2 c11 c12 x11 x12 c21 c22 x21 x22 … … c1j x1j c2j x2j … … c1n x1n c2n x2n a1 a2         Di  Dm ci1 ci2 xi1 xi2  cm1 cm2 xm1 xm2 …  … cij xij  cmj xmj …  … cin xin  cmn xmn ai  am bj b1 b2 … bi … bn ai bj Traducerea algebrică a celor trei tipuri de restricţii conduce la: (min)f(x) im1 jn1cijxij, x= x11, 12, jn1 xij ai, i 1, m im1 xij bj, j 1, n xij 0, i 1,m; j 1, n T mn x (6.1) Observaţie. Dacă disponibilul este mai mic decăt cererea sau mai mare, pentru a echilibra problema se introduce un depozit fictiv sau un centru fictiv cu cantitatea (fictivă) existentă sau cerută astfel încăt problema să devină echilibrată. Costurile unitare de transport pentru depozitul fictiv sau centrul fictiv se consideră nule. 66 j i S S S Teorema 6.1. Problema echilibrată de transport echilibrată (6.1.) admite cel puţin o soluţie admisibilă, iar o soluţie de bază admisibilă de bază are cel mult mn1 componente diferite de zero. Demonstraţie: Cum m a  n b S, considerând x aibj avem: i1 j1 i j   ij S n j 1 ij x m x i 1 ij n aibj = ai n b = ai S=a , i 1,m j=1 j=1 m aibi = bj m a = bj S=b , j 1,n i=1 i=1 i j S S S şi cum xij  0, ajungem la concluzia că X(xj ) verifică sistemul de restricţii şi condiţiile i nenegativitate, fiind deci soluţie admisibilă a problemei de transport. Putem observa că din cele mn restricţii doar mn1 sunt independente. Este suficient pentru aceasta, de exemplu, ca din suma primelor m restricţii, să scădem suma următoarelor n1 şi vom obţine atunci ultima restricţie. Pe de altă parte ştim că numărul componentelor nenule ale unei soluţii de bază admisibile de bază este cel mult egal cu numărul restricţiilor independente şi în felul acesta teorema este demonstrată. 6.1. SOLUŢII ADMISIBILE DE BAZĂ PENTRU O PROBLEMĂ DE TRANSPORT Soluţiile optime pentru o problemă de transport vor fi căutate în mulţimea soluţiilor admisibile de bază şi din acest motiv vom pune în evidenţă două metode prin care putem găsi aceste soluţii. a) Metoda diagonalei. Se începe cu determinarea componentei corespunzătoare căsuţei din stânga sus a tabelului deci 11 x Vom lua 11 min(a ,b ) şi vom atribui valoarea zero tuturor variabilelor de pe aceiaşi linie sau x 1 1 coloană cu 11 în funcţie de următoarele situaţii: x 1. dacă a1  b1, vom lua 11 a , 1j 0, j 2,n x x 1 67 2. dacă a1  b1 , vom lua 11 1, i1 0, i  2,m x x b 3. dacă a1  b1 , vom lua 11 a 1, iar 12 sau x21 0, x 1 b x şi restul componentelor de pe linia şi coloana lui 11 egale cu zero. În ultima situaţie obţinem x soluţii degenerate. Modificăm apoi a1 şi b1 înlocuindu-i cu 1 11 x 1, b 11 11. Procedeul se repetă a a 1 1 b x la pasul următor pentru calculele rămase necompletate care formează un tablou cu m – 1 linii şi n coloane ; m linii şi n 1 coloane ; sau m 1 linii şi n 1 coloane după cum primul pas s-a desfăşurat în situaţiile 1), 2), respectiv 3). Prima componentă ce va fi determinată acum va ocupa aceeaşi poziţie (colţul stânga – sus), în noul tabel. Procesul se termină în cel mult m n 1 paşi, la fiecare pas determinându-se complet o linie (situaţia 1), o coloană (situaţia 2), sau o linie şi o coloană situaţia 3. La fiecare pas se determină o singură componentă nenulă. Componentele nebazice nu se completează în tabel, căsuţele corespunzătoare rămânând libere pentru a nu se confunda componentele nebazice cu eventualele componente bazice nule. b) Metoda costului minim. Această metodă este întru-totul similară celei a diagonalei, cu singura diferenţă că se va determina mai întâi componenta xij pentru care cij minckl , luând . xij min(ai,bj) . Exemplul 6.1. Să se determine o soluţie admisibilă de bază pentru problema de transport D C C1 C2 C3 C4 ai D1 4 2 2 3 35 D2 2 2 4 5 40 D3 3 4 2 3 20 bj 50 10 25 10 95 95 68 a) Metoda diagonalei. Problema este echilibrată. Vom căuta o soluţie cu m n 1 3 4 1 6 35 4 2 2 3 35, 0 2 2 4 5 40,25, 15, 15 10 15 5 0 3 4 10 2 10 3 20, 10, 0 50 10 25 10 15 0 10 0 0 0 Pasul întâi: deoarece min(a1,b1)=a=35, punem conform situaţiei 1) 11 35, 12 13 14 0. Modificăm pe a şi b înlocuindu-i cu a a 11 0, x x x x 1 1 1 1 x b b 11 15. 1 1 x Pasul al doilea: se pleacă de la tabelul format prin înlăturarea primei linii,: deoarece min(a2,b1)=b1=15, punem (conform 2)) x21 15, x31 0 (x31 componentă nebazică). Modificăm pe a2 şi b înlocuind cu a2 a2 x21 25, 1 b x21 0. 1 b 1 Pasul al treilea: a rămas de completat tabelul format din liniile a doua, a treia şi ultimile trei coloane. Cum min(a2,b2)=b2=10 punem x22 10, x32 0 (x32 componentă nebazică) înlocuindu-i cu a2 a2 x22 15, b2 b2 x22 0. Pasul al patrulea: a rămas de completat tabelul format din ultimile două linii şi coloane. Cum min(a2,b3)=a =15, punem (conform 1)) x23 15, x4 0(x24 componentă nebazică). Modificăm pe a23 şi b3 înlocuindu-i cu a2 a2 x23 0, b b x23 10. 3 3 Pasul al cincilea: a rămas de completat tabelul format din ultima linie şi ultimile două coloane. Cum min(a3,b3)=10 punem x33 10. Modificăm pe a3 şi b3 înlocuindu-i cu a3 a3 x33 10, 3 b x33 0. b 3 69 12 13 1 Pasul al şaselea: obligatoriu vom pune x34 10. S-a obţinut o soluţie admisibilă de bază: 11 x x 33 35,x 0,x 10,x34 10 0,x 4 0,x21 15,x22 10,x23 15,x24 0,x31 0,x32 0, costul total al transportului fiind 300 u.m. b) Metoda costului minim. Vom considera din nou tabelul simplificat: 4 2 25 2 10 3 35, 10, 0 30 2 10 2 4 5 40, 30, 0, 20 3 4 2 3 20, 0 50 10 25 10 0 0 0 0 20 Pasul întâi: mincij=c12 c13 c21 c22 c33 2. Putem determina ori care dintre componentele 12, 13,x21,x22,x33 . Luăm x22=min(a2,b2)=10, 12 x32 0 (componente x x x nebazice). Înlocuim pe a2 şi b2 cu a2 a2 x22 30, b b x22 0. 3 2 Pasul al doilea: min cij= c 3= c21= c33= 2. Luăm x21=min(a2,b1)=30, x23 x24 0 1 (componente nebazice). Înlocuim pe a2 şi b2 cu a2 a2 x21 0, b b x21 20. 1 1 Pasul al treilea: min cij= c13= c33= 2. Luăm 13 min(a ,b ) 25, x33 0 (componentă nebazică). Înlocuim a1 şi b3 cu x 1 3 a a 13 10, b b 13 0. 1 1 x 3 3 x Pasul al patrulea: min cij= c14= c31= c33= 3. Luăm x31 min(a3,b ) 20, x 1 0 (componentă bazică), x34 0 (componentă 1 1 nebazică). Înlocuim a3 şi b1 cu a3 a3 x31 0, b b x31 0. 1 1 Pasul al cincilea: obligatoriu 14 10. Am obţinut pentru problemă o soluţie x admisibilă de bază degenerată ( 11 0 componentă bazică). Ea este indicată în tabelul de mai x sus.Costul unui asemenea transport este 220 u.m. Metoda costului minim ne conduce de cele mai 70 multe ori la o soluţie admisibilă de bază mai bună decât metoda diagonalei, în sensul că realizează o valoare a cheltuielilor de transport mai mică. 6.2. SOLUŢII OPTIME PENTRU O PROBLEMĂ DE TRANSPORT Soluţiile optime (soluţii admisibile ce realizează minimul funcţiei obiectiv) ale unei probleme de transport le vom căuta în mulţimea soluţiilor admisibile de bază nedegenerate. Să considerăm că am determinat pentru problemă o soluţie admisibilă de bază nedegenerată: x=(xj)T . Pornind de la aceasta să vedem cum putem obţine o soluţie admisibilă. i Convenim ca o pereche de indici (i, j) să o numim celulă. Când vom spune celulă liberă înţelegem că ea se referă la o componentă nebazică (secundară), adică nu se găseşte printre cele m+n-1. Celule ocupate vor fi cele care se referă la componentele bazice (principale). Să considerăm o celulă liberă (i, j). Definiţia 6.1: Un şir de celule care începe cu celula liberă (i, j) şi se termină cu aceasta, conţinând în rest numai celule ocupate şi anume câte două din aceiaşi linie şi coloană se numeşte ciclul celulei libere (i, j). Cel mai simplu ciclu al unei celule libere (i, j) se prezintă sub forma de mai jos unde în dreptul celulelor am scris şi componentele corespunzătoare ale soluţiei. xij (i, j) (i1j ) 1 1 (i1,j) (i1,j ) 1 Să atribuim componentei xij (componenta nebazică), o valoare 0.Vom obţine atunci 1 i j x 1 1 i j x un vector care are mn componente diferite de zero, acesta trebuind să verifice sistemul de restricţii. 71 1 1 x x 1 i j i Acestea sunt verificate dacă vom scădea şi aduna succesiv valoarea  la componentele care au intrat în ciclul celulei libere (i, j). Indicăm pe acelaşi ciclul acest lucru. ij x 1 (i, j) (i1j ) 1 (i1,j) (i1,j ) 1 i j i j 1 x 1 x 1 Pentru a fi verificate condiţiile de nenegativitate va trebui ca să fie luat astfel încât ij1 0 şi x1 j 0. x i Odată îndeplinite aceste condiţii putem spune că am obţinut o nouă soluţie admisibilă x(xj ) unde: i xij , xij ij xi j i j , xke , x1 1 x1 j xke pentru k i,i ; j j, j 1 1 1 (6.2) Dorind acuma ca noua soluţie admisibilă să fie şi admisibilă de bază este suficient să considerăm: min(xj ,x1 j ) (6.3) 1 i i În acest fel se determină care componentă bazică devine nebazică. Relaţia (6.3) reprezintă forma criteriului de ieşire din algoritmul pentru găsirea soluţiilor optime a unei probleme de transport. Să vedem acum care sunt condiţiile ca soluţia admisibilă de bază x să fie soluţie optimă. Să evaluăm diferenţa f x -f x . Vom avea: f x -f x sau m n k 1e 1 ke c xke xke ij ij i j i j c c 1 1 1 c 1 c f x f x ij (6.4) unde 72 ij  cij cij ci1j ci1j (6.5) 1 1 care reprezintă suma algebrică cu semn alternat a costurilor celulelor corespunzătoare ce au intrat în ciclu şi care începe cu cij . Din (6.5) se observă că dacă pentru toate celulele libere (i, j) vom avea ij  0, şi aşa cum 0, soluţia admisibilă de bază x(xj ) considerată reprezintă soluţia optimă a problemei de i transport. Rezultă în felul acesta că criteriul de optimalitate este ij  0 (6.6) pentru toate celule libere (i, j). Tot din (6.4) se constată că dacă există δij aceasta fie: δij=maxδke ke δ 0 0, soluţia x poate fi îmbunătăţită. Pentru (6.7) Atunci componentei nebazice xij îi vom atribui o valoare 0. Relaţia (6.7) reprezintă forma criteriului de intrare în algoritm. Mărimile δij joacă rolul diferenţelor zj cj din problemele de programare de care ne-am ocupat. Sintetizând cele de mai sus în practică procedăm în felul următor: a) Se determină o soluţie admisibilă de bază nedegenerată prin una din cele două metode prezentate. b) Pentru toate celulele libere se calculează cantităţile ij . - dacă toate cantităţile ij sunt mai mici sau egale cu zero, atunci soluţia este optimă; - dacă există ij  0 se trece la etapa următoare. c) Se aplică criteriul de intrare. Vom atribui o valoare 0 componentei xij (nebazică) cea pentru care avem îndeplinită relaţia (6.7). d) Se consideră ciclul celulei xij , se introduce în ciclu şi se aplică criteriul de ieşire. Se determină valoarea lui după relaţia (6.3) indicându-se deci care componentă bazică devine nebazică. e) Etapele algoritmului se reiau până se găseşte soluţia optimă a problemei. 73 Observaţii: a) Prezenţa semnului de egal în criteriul de optim şi de intrare în bază, indică existenţa mai multor soluţii optime ; b) Pe parcursul aplicării algoritmului pot să apară soluţii admisibile degenerate care pot conduce la un fenomen numit fenomen de ciclaj. Pentru evitarea apariţiei soluţiilor admisibile degenerate se întrebuinţează aşa-numita metodă a perturbării. Aceasta constă în înlocuirea cantităţilor ai cu ai , şi cantităţii bn cu b m (unde reprezintă o cantitate foarte mică, strict pozitivă). La finalul algoritmului, după determinarea soluţiei optime a problemei peturbate, se face  0, obţinându-se soluţia optimă a problemei iniţiale. n c) Dacă la aplicarea criteriului de intrare sunt mai multe mărimi δij pozitive maxime egale putem lua oricare dintre acestea. Exemplul 6.2. Să se determine soluţiile optime ale problemei de transport din exemplul 6 Pentru a fi siguri că nu apar soluţii admisibile de bază degenerate aplicăm de la început metoda perturbării, adică considerăm: a1=35+ε, a2=40+ε, a3=20+ε, b1 50, b2 10 b3 25, b4 10, 3ε O soluţie admisibilă de bază determinată prin metoda costului minim din tabel este dată în tabelul: D C C1 C2 C3 C4 4 2 2 3 1 10 25 ε D 2 2 3 5 2 40+ε D 3 4 2 3 3 10 ε 10+2ε D ai 35+ε 40+ε 20+ε bj 50 10 25 10+3ε 95+3ε 95+3ε Să găsim cantităţile ij pentru celulele libere. Considerăm ciclul celulei libere (1,1) şi în dreptul celulelor scriem costurile. 74 4 3 (1,1) (1,4) (3,1) (3,4) 3 3 Vom avea: δ11= 4 3 3 3 1 Analog: δ22 23 24 32 33 δ δ δ δ 2 2 3 3 3 2 1 4 2 3 3 3 2 3 5 3 3 2 3 4 2 3 3 2 2 2 3 3 0 Deoarece toate cantităţile ij sunt mai mici sau egale cu zero am găsit soluţia optimă a problemei. Făcând ε=0 găsim: 11 0, x 2 10, x 25, x 4 0, x21 40, x22 0, x23 0, x24 0 x 1 1 1 x31 10, x32 0, x33 0, x34 0, (min)f(x) 210 u.b. Deoarece δ33 0 , algoritmul poate continua atribuind componentei x33 (nebazică), o valoare 0. Considerăm ciclul celulei (3,3) şi introducem pe în ciclu pentru a satisface restricţiile. 25 (1,3) (3,3) Aplicăm criteriul de ieşire luând ε (1,4) (3,4) min(25, 10 2ε)=10+2ε . 10 2ε 75 Componentele noii soluţii sunt date în tabelul: D c C1 C2 C3 C4 4 2 2 3 1 10 15-2ε 10+3ε D 2 2 4 5 2 40+ε D 3 4 2 3 3 10 ε 10+2ε D ai 35+ε 40+ε 20+ε bj 50 10 25 10+3ε 95+3ε 95+3ε Făcând ε=0 găsim soluţia optimă. 11 0, x 2 10, x 3 15, x 4 10, x21 40, x22 0, x23 0, x24 0 x 1 1 1 x31 10, x32 0, x33 10, x34 0, (min)f(x) 210 u.b. Exemplul 6.3 Să se determine soluţiile optime ale problemei de transport: c D D1 D2 C1 C2 C3 ai 1 2 3 50 3 4 1 25 bj 40 20 30 90 75 Problema nu este echilibrată. Vom adăuga un depozit cu costuri zero. Vom avea deci: c D D1 D2 D3 C1 C2 C3 ai 1 2 3 50 3 4 1 25 0 0 0 15 bj 40 10 25 90 90 76 O soluţie admisibilă de bază determinată prin metoda costului minim din tabel este prezentată în: D c C1 C2 C3 ai D1 25 1 20 2 5 3 50 D2 3 4 25 1 25 D3 15 0 0 0 15 bj 40 20 30 90 90 Componentele soluţiei au fost determinate în ordinea: x31 15, x 1 25, x21 25, x 2 20, 1 1 13 5, x21 0, x22 0, x32 0, x33 0 x Pentru celulele libere calculăm cantităţile ij . δ21 22 32 33 δ δ δ 3 1 3 1 4 4 2 3 1 4 0 2 1 0 1 0 3 1 0 2 Soluţia nu este optimă. Avem că 33 max ij 2. Considerăm ciclul celulei (3,3) ij 0 componentei x33 (nebazică) îi atribuim o valoare 0 şi introducând pe în circuit avem: 25 (30) 5 (0) (1,1) (1,3) (3,1) (3,3) 15 (10) (5) După criteriul de ieşire luăm min(1,15) 5. Noua soluţie este prezentată în tabelul: 77 D c C1 C2 C3 ai D1 30 1 20 2 3 50 D2 3 4 25 1 25 0 0 0 3 10 5 D 15 bj 40 20 30 90 90 Pentru celulele libere vom avea: δ13 21 22 32 δ δ δ 3 1 0 0 2 3 1 0 0 2 4 2 1 0 0 1 2 0 2 1 0 1 Soluţia nu este optimă. Atribuim componentei x32 o valoare 0. 30 (40) 20 (10) (1,1) (1,2) (3,1) (3,2) 10 (0) (10) Noua soluţie este prezentată în tabelul: D c C1 C2 C3 ai D1 40 1 10 2 3 50 D2 3 4 25 1 25 0 0 0 3 10 5 D 15 bj 40 20 30 90 90 78 Această soluţie este optimă deoarece avem: 13 1, 21 3, 22 3, 31 1. Soluţia problemei iniţiale va fi: 11 40, 12 10, x 3 0, x21 0, x22 0, x23 25.(min)f(x) = 85u.b. x x 1 Se observă că cerinţele centrelor C2 şi C3 nu vor fi satisfăcute în totalitate. 79 US7 Serii de numere reale Problema sumelor infinite (a seriilor) este adesea întâlnită în asigurări și finanțe. Rezolvarea ei presupune a stabili dacă o astfel de sumă există și în caz afirmativ, care este valoarea ei. Nu este o problemă ușor de rezolvat pentru ca nu toți studenții cunosc bazele analizei matematice (mai precis teoria șirurilor de numere reale). Prin activitățile tutoriale completăm acest lucru, după care definim noțiunile fundamentale (serie convergentă, serie divergentă, suma unei serii convergente, criteriu de convergență etc.). Se insistă asupra faptului că suma unei serii convergente se poate aproxima, după ce (printr-un criteriu dintr-o listă destul de bogată) am stabilit natura ei. Studenții trebuie să înțeleagă mai întâi ce înseamnă serie convergentă analizând exemplele clasice: seria geometrică, armonică, armonică generaliată, telescopică. Urmează apoi reținerea criteriilor care se aplică oricărei serii (criteriul general de convergență, criteriul general de divergență). După această etapă urmează criteriile pentru fiecare tip de serie (criterii pentru serii cu termeni pozitivi, criterii pentru serii cu termeni alternați, criterii pentru serii cu termeni oarecare). Deși este vorba de un singur subiect teoretic, sunt necesare aproximativ 8 ore pentru însușirea lui: două ore pentru noțiunile teoretice generale, două ore pentru seriile remarcabile, două ore pentru seriile cu termeni pozitivi și două ore pentru seriile cu termeni oarecare și pentru seriile alternate. Pozițiile [1] și [3]-[6] din lista bibliografică pot fi de un real ajutor, dar pot fi utile si alte tratate de analiză matematică. * * * Este bine cunoscut faptul că, în prezent, fenomenele economice devin din ce în ce mai complexe şi interdependente; mai nou, tendinţa de globalizare a economiei (şi nu numai) duce la o concurenţă acerbă pe mai toate pieţele şi în mai toate domeniile de activitate economică, fenomen care se simte din ce în ce mai pregnant şi în economia românească. Aceste situaţii, cer din partea specialiştilor în economie, cunoştinţe cât mai precise în vederea observării, înţelegerii, modelării şi rezolvării pe baze riguros ştiinţifice a problemelor 80 complexe cu care se confruntă. Această abordare nu mai este astăzi posibilă, fără recurgerea la un aparat matematic cât mai adecvat. Modul mai puţin riguros, chiar empiric, de luare a deciziilor economice doar pe baza unei simple experienţe, este total depăşit, astfel încât cei care nu se vor adapta noilor condiţii şi cerinţe economice, vor fi eliminaţi încet, încet de pe piaţă. În acest capitol se prezintă, într-o manieră accesibilă studenţilor economişti, câteva elemente de bază din teoria seriilor numerice şi a seriilor de funcţii (cazul particular al seriilor de puteri). Aplicaţiile practice ale acestei teorii se întâlnesc în multe domenii de cercetare din matematică, fizică, chimie, biologie etc. În particular, aplicaţiile cu caracter economic ale teoriei seriilor sunt multiple şi permit modelarea şi rezolvarea unor probleme economice cu caracter atât determinist cât şi aleator (probabilistic). Elemente şi aplicaţii ale teoriei seriilor sunt regăsite în Teoria probabilităţilor şi Statistica economică, Micro- şi Macroeconomie, Gestiunea stocurilor şi a fenomenelor de aşteptare etc. Pentru parcurgerea şi înţelegerea noţiunilor prezentate în acest capitol sunt necesare noţiuni de teoria limitelor şi a şirurilor studiate la Analiza Matematică din clasa a XI-a. 81 n ; 7.1 SERII NUMERICE 7.1.1. ŞIRURI DE NUMERE REALE Definiţia 7.1.1: Numim şir de numere reale, o aplicaţie: f : f(n) an Vom nota şirul cu simbolul: (an)n ;an se numeşte termenul general al şirului. Definiţia 7.1.2: Spunem că şirul (an)n este convergent către un număr a , dacă, orice vecinătate a lui a conţine toţi termenii şirului cu excepţia unui număr finit de termeni. Vom nota cu an a sau lima a . n Numărul "a" se numeşte limita şirului şi dacă există este unică. Dacă şirul (an)n nu este convergent, atunci acesta se numeşte şir divergent. Un şir divergent sau nu are limită sau limita sa este sau . Exemple: a) an n n2 1 bn n2 2n2 3 ; cn 1 1 n n sunt şiruri convergente (cu limitele 0, 1 2 respectiv e=2,7.11…; b) an n ; bn n3n 1 ; cn 3 1 n2 sunt şiruri divergente cu limitele , , respectiv ; 2 c) an ( 1)n, bn ( 1)n ln(n2 n 1) ,cn sin(n 1) sunt şiruri divergente care nu au limită. 82 Teorema 7.1.1 (de caracterizare a şirurilor convergente) Şirul an a dacă şi numai dacă: ( ) 0,( )n astfel încât an a , pentru ( )n n . Teorema 7.1.2. (Cauchy) Şirul (an)n este convergent dacă şi numai dacă: ( ) 0, ( )n astfel încât : an k an , pentru ( )n n şi ( )k * . Observaţii: i) Cele două teoreme, ca şi definiţia, nu sunt folosite în mod curent în aplicaţii practice; cu ajutorul lor se pot demonstra rezultate a căror aplicare este mult mai simplă (de exemplu, teorema de convergenţă a şirurilor monotone şi mărginite, teorema cleştelui etc.). ii) Aşa cum se observă, aplicarea teoremei 7.1.2. nu presupune cunoaşterea apriorică a limitei a, în enunţul acesteia intervenind numai termenul general an respectiv an k ; acest lucru este deosebit de util în aplicaţiile care conţin şiruri recurente sau în care limita şirului este dificil de determinat. 7.1.2 NOŢIUNEA DE SERIE NUMERICĂ. CONVERGENŢĂ. PROPRIETĂŢI GENERALE ALE SERIILOR. După cum se ştie, oricărei mulţimi finite a1,a2,,an de numere reale i se poate asocia un număr real s  a1 a2 an numit suma acestora. Se pune întrebarea firească: dacă  mulţimea conţine o infinitate de numere reale a1,a2, ,an, (cu alte cuvinte un şir de numere reale (an)n ) i-am putea asocia un număr S care să generalizeze conceptul de sumă finită? Fie şirul : (an)n :a0,a1,a2, ,an, 83 1 Vom asocia acestuia, un nou şir de numere reale, (Sn)n astfel: S0  a0   S1  a0 a1  S2  a0 a1 a2   .......................... Sn  a0 a1 a2 an  k ak  n 0   :S0,S ,S2, ,Sn, definit (7.1.1) pe care îl vom numi şirul sumelor parţiale ataşat şirului (an )nN . Definiţia 7.1.3. Numim serie numerică, de termen general an , expresia n an  lim Sn  a0 a1 a2 an  (7.1.2) def   0 n Observaţie: În unele cazuri numerotarea termenilor seriei poate începe de la n 1  a  sau n1       n de la n  n oarecare fixat   a  (vezi şi proprietăţile seriilor numerice).    n 0 nn0  Definiţia 7.1.4.: Spunem că seria a este o serie: n 0 n a) convergentă (C ), dacă şirul (Sn)n b) divergentă (D) , dacă şirul (Sn)n este şir convergent; este şir divergent. Cu alte cuvinte problema convergenţei unei serii numerice este echivalentă cu problema convergenţei unui şir de numere reale şi anume, şirul sumelor parţiale (Sn)n . Observaţii: i) Dacă S S, vom spune că seria a (C) şi are suma S (notăm: a S); n 0 n 0 n n n ii) Dacă S , spunem că seria a (D) şi are suma  notăm a ); n 0 n 0 n n n 84 n iii) Dacă (Sn)n nu are limită, spunem că a (D) şi nu are sumă. n 0 În studiul seriilor numerice, trebuie cercetat (ca şi în cazul şirurilor) două aspecte: - natura seriei (dacă este convergentă sau divergentă) - suma seriei (în caz de convergenţă) Dacă pentru şiruri stabilirea convergenţei şi determinarea limitei acestora este în cele mai multe cazuri una şi aceeaşi chestiune, în cazul seriilor situaţia se schimbă radical; de cele mai multe ori este relativ uşor să stabileşti convergenţa seriei, dar foarte greu, ba câte odată chiar imposibil de determinat suma acesteia. De obicei, cu metodele specifice analizei numerice, se determină o aproximaţie S a sumei S a seriei (S S). Pentru a se cerceta natura unei serii, nu se poate folosi întotdeauna definiţia, ci adesea se utilizează condiţii echivalente cu aceasta sau numai condiţii suficiente de convergenţă numite “criterii de convergenţă”. Vom prezenta unul dintre aceste criterii în cele ce urmează. Teorema 7.1.3. (criteriul general de convergenţă a lui Cauchy) Seria a este (C) dacă şi numai dacă: ( ) 0,( )n astfel încât: n 1 n an 1 an 2 an k , pentru ( )n n , ( )k . Această teoremă se obţine ca o consecinţă a Teoremei 7.1.2, considerând în loc de şirul (an)n , şirul sumelor parţiale (Sn)n . Într-adevăr, conform (7.1.1) avem: Sn k Sn (a1 an an 1 an k ) (a1 an) an 1 an k 85 ; 7.1.3. EXEMPLE REMARCABILE DE SERII 1) Seria geometrică aqn; cu a,q n 0 Termenul general al acestei serii are forma an aqn (*) şi depinde de doi parametri (constante reale) a, q. Evident, este logic ca natura seriei şi suma acesteia în caz de convergenţă să depindă de aceşti parametri. Astfel, dacă: i) a 0, atunci an 0, ( ) q , deci seria este convergentă şi are suma S 0. ii) a 0, avem: (*) Sn a0 a1 a2 an a aq a(1 q q2 qn) aq2 aqn deci Sn reprezintă suma unei progresii aritmetice cu raţia q (de aici şi numele seriei) şi conform rezultatelor de la acest capitol (vezi manualul de Algebră, clasa a X-a) vom avea: 1 qn 1 a ; Sn 1 q q 1 a(n+1); q=1 Trecând la limită în relaţia de mai sus, obţinem: a 1 q q ( 1,1) lim S ; q [1, ) (depinde de semnul lui a) n n ( ) ; q ( , 1] cu alte cuvinte, pentru: q 1,1 aqn n 0 q ( ,1] [1, ) 1 aq seria converge cu suma S 1 aq (7.1.3) aqn(D) cu suma S sau nu există n 0 86 Observaţie: De fapt, această serie cuprinde o infinitate de serii particulare (pentru a, q având valori fixate) a căror natură o cunoaştem; foarte important este că, în caz de convergenţă, ştim şi suma seriei. Exemple: a) 1 n 0 2n n , este serie geometrică cu a 1, q n 0 1 2 1 şi conform (7.1.3) rezultă că: 2 1 n 0 2n 1 1 1 1 2 22 2n 2 (C) b) ( n 0 1)n 9n n 0 n 3 , este serie geometrică cu a 1 şi q 1 1 3 şi conform (7.1.3) rezultă că: ( 1)n n 0 9n 1 1 1 ( 1)n 3 9 81 9n 4 (C) c) 5n 1 n 0 22n 5 n 0 n 4 , este serie geometrică cu a 5, q 5 5 4 şi conform (7.1.3) seria n 1 5 diverge în acest caz n 0 22n . d) ( 1)n 1 2 n ( 1) 2 n este serie geometrică cu a 1, q= 2 şi deci n 0 n 0 seria este divergentă; mai mult suma nu există. 2) Seria armonică generalizată 1 ; R n 1 n Evident, ca şi în şirul seriei geometrice, natura (şi în caz de convergenţă suma) seriei armonice generalizate, va depinde de valoarea parametrului . Aplicând, Teorema 4.2 se poate arăta că, seria: (C) ; 1 1 n 1 n (D) ; 1 (7.1.4) 87 Observaţii: i) În caz de convergenţă ( 1), nu există o formulă generală de stabilire a sumei S în funcţie de parametrul , asa cum s-a stabilit în cazul seriei geometrice; ii) Evident, în cazul în care seria diverge ( 1), suma seriei este S (ca serie cu termeni pozitivi). Exemple: a) n 1 1 n (D),deoarece 1 n 1 1 n 1 1 1 1 2 3 n ) b) n 1 1 n2 (C),doarece 2 1 n 1 1 n2 1 1 1 1 22 32 n2 S c) n 1 1 1 3 n2 n 1 n2 3 (D), deoarece 2 3 1 Seria n 1 1 n se numeşte seria armonică (simplă) sau seria lui Riemann, deoarece termenul general an este media armonică a vecinilor săi, adică satisface relaţia 2  1  1 . n n1 n1 a a a 3) Seria telescopică n 1 1 n(n 1) Această serie este convergentă. Într-adevăr scriind termenul general an  n(n 1) sub 1 forma (*) an 1 n n1 1, obţinem termenul general al şirului sumelor parţiale (Sn)nN de forma: Sn  a1 a2 an  12  23  34  n(n 1)  1 1 1 1 (*) 1 1 1 1 1 1 1 2 2 3 3 4 n n1 1 = (după simplificări) = 1 n 1. 1 88 , Atunci n 1 1 n(n 1) lim S lim 1 n n n 1 n 1 deci n-1 1 n(n 1) 1 (C) Dacă notăm cu n 1 , atunci termenul general al seriei an  n(n 1) , se poate scrie 1 n sub forma: an  n n1 (7.1.5) Seriile numerice a căror termen general an au forma de mai sus, se numesc serii telescopice şi se poate observa uşor ca seria n an (C)  (n)nN (C).   1 Mai mult seria va avea suma S 1 , unde 1 este primul termen al şirului (n)nN , iar   lim  . n n Vom enunţa în continuare câteva din proprietăţile generale ale seriilor numerice, care se pot demonstra relativ uşor, ţinând cont de definiţia seriei şi a operaţiilor cu şiruri convergente. 7.1.4. PROPRIETĂŢI GENERALE ALE SERIILOR NUMERICE (P1) Dacă unei serii  a , i se adaugă sau i se suprimă un număr finit de termeni, n1  n natura acesteia nu se schimbă. Observaţie. Evident, în caz de convergenţă, suma S a unei serii se modifică corespunzător prin adăugarea (sau scăderea) a sumei termenilor finiţi care se adaugă (sau se suprimă). (P2) Sacă unei serii convergente  a , i se asociază termeni în grupe finite, cu n1  n păstrarea ordinii lor, atunci se obţine tot o serie convergentă şi cu aceeaşi sumă. 89  n  n  n  n Observaţii: i) Proprietatea de mai sus, poate fi scrisă şi sub forma: (P ) Dacă n an  a1 a2 an S, atunci şi: 1   2 (a1 a2 an1 ) (an1 1 an2 ) (an2 1 an3 ) S ii) Asociind termenii unei serii divergente, în grupe finite cu păstrarea ordinii se pot obţine serii convergente. Exemplu: Fie seria ( 1)n 1 1 1 1 n 0 Atunci: ( 1)n (termenul general fiind an ( 1)n ) S0 1, S1 11 0, S2 1111, S3 1111 0, deci: S2n 1 şi S2n1  0 ; evident şirul sumelor parţiale (Sn )nN diverge deoarece are două subşiruri care converg la limite diferite. Rezultă că seria  (1)n (D) şi nu are sumă. n0  Dacă asociem termenii în grupe de câte doi (sau de câte patru, şase, etc.) obţinem: (11)(11)(11)(11) 000 0 deci noua serie obţinută converge şi are suma S = 0. (P3) Dacă seriile  a n1 şi  b sunt convergente cu sumele S1 respectiv S2, atunci seria n1  (a  b ) converge şi are suma S1 + S2. n1  n n Cu alte cuvinte dacă: n1 an S1,n1 bn  S2  n1 (an bn )  S1 S2 (7.1.6)       (P4) Dacă * , atunci seriile  a n1 şi  a au aceeaşi natură; mai mult în caz de n1 convergenţă, dacă  a  S atunci  a  S. n1 n1 n n   90 n n n Teorema 7.1.4. Dacă seria a (C), atunci lim a 0 n 1 n n n Demonstraţie: Fie Sn  a1 a2 an ,()nN*. Deoarece seria converge, avem: n1 an  lim Sn  S (1)   n  Dar an Sn Sn1; trecând la limită în această relaţie obţinem: (1) lim an  lim(Sn Sn1)  lim Sn  lim Sn1 SS  0 n n n n q.e.d. Observaţie: Condiţia din teorema de mai sus (an 0) este o condiţie necesară dar nu şi suficientă pentru convergenţa unei serii; deci, dacă termenul general an 0 acest lucru nu implică obligatoriu că seria  a (C). Ca exemple se pot da cazurile particulare ale seriei n1  n armonice. Astfel  1 (D) chiar dacă a  1  0. n1  n n n O formulare echivalentă a teoremei 7.1.4. des folosită în aplicaţii practice, este: Teorema 7.1.5. Criteriul general de divergenţă (C.G.D) Fie seria a , n 1 an . Dacă lim a n 0 atunci seria a (D). n 1 Exemple: a) n1 2n 1 (D) , deoarece lim an  lim 2n 1  2  0 n n 1   n n b)  n 1n (D), deoarece: n1    n   lim an  lim n 1n  lim 1 1n  e  2,71 0 n n         n n n c)  1 n n2  n1 2n (D) , deoarece lim a lim 1 n n n n n2 2n 0 91 n n  n n n În concluzie, dacă: i) lima 0 n an (D) n 1 ii) lima 0 n acest caz) an poate fi (C) sau (D) (trebuie făcută o analiză suplimentară în n 1 7.2 SERII CU TERMENI POZITIVI Definiţia 7.2.1. Numim serie cu termeni pozitivi , o serie de forma:  a cu a 0, ( )n . n1  n n Teorema 7.2.1. (de caracterizare a convergenţei seriilor cu termeni pozitivi) Seria cu termeni pozitivi  a , a  0 converge dacă şi numai dacă şirul sumelor n1  n n parţiale (Sn)n este mărginit. Demonstraţie: Avem Sn1 Sn  (a1 a2 an1)(a1 a2 an1) an1  0deci Sn1 Sn , adică şirul (Sn )nN este monoton crescător. Dacă (Sn)n este şi marginit (superior), rezultă că este convergent (conform teoremei de convergenţă aşirurilormonotoneşi mărginite)deci  a (C). n1  n Reciproc, dacă  a (C)  (S ) n1 deci teorema este demonstrată. (C). Dar un şir convergent este obligatoriu mărginit, q.e.d. 92 n n n Observaţie: În aplicaţiile practice, Teorema 7.2.1. nu este folosită foarte des deoarece, în general, este relativ complicat a se demonstra mărginirea şirului (Sn)n ; întradevăr, trebuie demonstrată relaţia: Sn a1 a2 an M, ( )n * cu M 0constantă care nu depinde de n . Mult mai utile sunt următoarele rezultate (numite criterii de convergenţă), care pot fi clasificate în două categorii: criterii de comparaţie şi criterii cu limită. I) Criteriul de comparaţie de specia I (cu inegalităţi simple) Fie seriile: a şi b care satisfac relaţiile 0 n 1 n 1 Atunci: an bn, ( )n . i) Dacă b (C) a (C) n 1 n 1 n n ii) Dacă a (D) b (D) n 1 n 1 n n Exemplu: Fie seria n! . Termenul general a n! satisface: n 1 n n n n n 0 an 1 2 3 n 1 2 3 n n n n n n n n 1 2 2 not n n n n2 n b Dar: b n 1 2 n 1 n2 2 1 (C) (ca serie armonică generalizată 2 n n 1 n 1 1 n cu 2 1) Atunci conform punctului i) din criteriu n! (C). n 1 n n 93 n n II) Criteriul de comparaţie de specia a III-a (cu limită) Fie seriile a şi b cu a , b 0. Dacă există: n 1 n 1 n n n n lim an k (0, ) (7.2.1) b n n atunci seriile au aceeaşi natură. Exemplu: Fie seria n 1 1 2n 3 cu an 1 2n 3 0 Luăm ca serie de comparaţie, seria armonică (simplă) b n 1 că diverge. Deoarece: 1 , despre care ştim n 1 1 lim an lim 2n 3 b 1 n n n n lim 2nn 3 n 1 2 k 1 2 (0, ) deci conform criteriului, seria a are aceeaşi natură cu seria b , deci diverge. n 1 n 1 n n Observaţii: i) Criteriul de comparaţie la limită, tratează şi cazurile când k=0 sau k , dar s-a renunţat la prezentarea lor pentru simplificarea expunerii. ii) Aplicarea criteriilor de comparaţie, nu este în general foarte simplă, deoarece trebuie satisfăcute simultan următoarele două condiţii: a) alegerea seriei de comparaţie (de exemplu, b ), funcţie de seria dată (de n 1 n exemplu, a a cărei convergenţă vrem să o stabilim); n 1 n b) trebuie cunoscută aprioric natura seriei b (convergentă sau divergentă n 1 n Prima condiţie este influenţată foarte mult de abilităţile matematice ale rezolvitorului, iar cea de a doua de “experienţa” sa. Pentru această din urmă condiţie, seriile geometrică şi 94 n ; armonică pot furniza o infinitate de serii particulare (a căror natură o cunoaştem) posibil de folosit ca serii de comparaţie. iii) În schimb, următoarele criterii (numite criterii cu limită) se aplică mult mai uşor, formularea acestora depinzând exclusiv de termenul general an al seriei an a n 1 cărei natură vrem să o determinăm. III. Criteriul raportului (D`Alembert) Fie seria a , a n 1 n 0 şi: lim an 1 0, . Atunci, dacă: n n a i) 1 a (C) n 1 n ii) 1 a (D) n 1 n iii) 1 ; caz de nedeterminare (nu se poate preciza natura seriei) Exemple: 1)  n! ; cu a  n! n1  n n n n n respectiv an1  ((n 1)!1 . Atunci n n 1)   lim a1  lim (n 1))!1  n!  lim (n 1)n(n)1)  n!  n n (n 1 n!( 1 n n n n n n n n a n (după simplificări) lim n n n 1 n 1 e 1, deci seria n! converge. n 1 n n 2) n 1 2n n2 1 n n1 cu an n2 1 şi an1  (n 1)2 1. 2 2 Atunci: lim n an 1 an lim n 2n 2 n2 1 n2 2n 2 2n 2 lim n n2 1 n2 2n 2 2 1 2 1deci, seria n 2 n 1 2   , diverge. n1 95  n n  n 2  3) n1 n ; cu an  n şi an1  n1 1. Avem: 1 1    lim an1  lim nn 1 1, deci conform acestui criteriu nu putem preciza natura seriei  a n n n  armonice  1 (despre care ştim însă că este divergentă, vezi formula (7.4)). n1  n IV. Criteriul rădăcinii (Cauchy) Fie seria  a , a  0 şi n1 lim n an   0,. Atunci, dacă:   n i) 1 a (C) n 1 n ii) 1 a (D) n iii) 1 ncazdenedeterminare(nusepoateprecizanaturaseriei) 1 Exemple: 1) n 1 n2 2n ; cu an n2 2n 0. Atunci: lim n a lim n n2 n n n n 2 lim (n n)2 2 n 1 2 2 1 deci, seria  converge. n1 n  2n 2) n 1 nn 1 n2 ; cu an n 1 n2 n 0. Atunci: lim n a lim n n n n n 1 n2 n lim n n 1 n n lim 1 n 1 n n e 1 deci seria diverge. 3)  1 n1 n2 ; cu an  1  0. Atunci:   lim n an  lim n n2  n  1 2  1 1, deci nu putem preciza cu acest criteriu, 1 lim  n 1 n n n natura seriei  1 (despre care ştim că este convergentă, ca seria armonică generalizată, n1  2 n  1 n1 n cu   2 1). 96 Observaţii: i) În unele din exemplele de mai sus am folosit limitele fundamentale: lim n n 1 şi n lim n P (n) 1 (7.2.2) n k unde: P (n) aknk ak 1nk 1 a2n2 a1n a0 polinom de grad k cu k coeficienţi reali şi ak 0 ii) Se observă că aplicarea acestor criterii, constă doar în calculul unei limite; totuşi dacă termenul general al seriei an , are o formă mai complicată determinarea lui poate fi o problemă deloc uşoară. Se poate apela în aceste cazuri la criteriile de comparaţie sau la alte criterii de o formă mai complicată care nu sunt prezentate aici (se poate consulta în acest scop, orice curs de analiză matematică care tratează convergenţa seriilor numerice). iii) Datorită relaţiei: lim an 1 λ lim n a (7.2.3) n n n a n rezultă că dacă un criteriu ne conduce la cazul de excepţie 1, nu se încearcă aplicarea celuilalt criteriu, deoarece se va obţine aceeaşi concluzie; obişnuit aplicarea criteriului raportului sau rădăcinii se face ţinând cont de forma termenului general an (astfel încât calculul limitei să fie cât mai simplu). iv) Se observă, că există situaţii – chiar pentru serii foarte simple – în care aceste criterii nu pot preciza natura seriilor (λ 1). De obicei în aceste cazuri se aplică criteriile de comparaţie respectiv criteriul general de divergenţă (C.G.D), sau următorul criteriu (a lui Raabe-Duhamel). v) V) Criteriul lui Raabe – Duhamel Fie seria a , a n 1 n n 0 şi lim n an 1 . Atunci, dacă: n n 1 a i) 1 a (C) n 1 n ii) 1 a (D) n 1 n iii) 1 ; cazdenedeterminare (nusepoateprecizanaturaseriei) 97 2 n Exemple: 1) n 1 1 n2 ; cu an 1 n2 0. Am văzut cum criteriul rădăcinii ne conduce la cazul de nedeterminare λ 1. Vom aplica criteriul lui Raabe – Duhamel pentru a “scăpa” din această situaţie. Avem: lim n an 1 n n 1 a lim n (nn 1)2 1 lim n (n2 n 2n 1) n2 n2 lim 2nn 1 2 1, deci seria converge. n 2) n 1 2n1 1 ; cu an 1 2n 1 0. Se verifică imediat că atât criteriul raportului cât şi cel al rădăcinii ne conduc la cazul λ 1. Dacă aplicăm criteriul lui Raabe – Duhamel obţinem: lim n an 1 lim n 2n n n 1 n a 2n 3 1 1 lim n 2n 2n 1 1 deci (în acest caz) nici criteriul lui Raabe – Duhamel nu ne “scoate” din cazul de nedeterminare. Dacă aplicăm în schimb criteriul de comparaţie la limită, considerând seria armonică : n 1 bn n 1 1 (D), obţinemk lim bn lim 2nn 1 2 . n a 1 n n n Deci cele două serii au aceeaşi natură, cu alte cuvinte seria n 1 Obeservaţii: 2n1 1 (D). i) Criteriul lui Raabe – Duhamel nu se aplică (în mod obişnuit) decât atunci când criteriul raportului sau rădăcinii ne conduc la cazul de nedeterminare λ 1. În cele mai multe dintre situaţii, acest criteriu ne scoate din cazul de nedeterminare; în caz contrar se va încerca aplicarea criteriilor de comparaţie. ii) Totuşi, există situaţii unde datorită formei termenului general al seriei aplicarea directă a criteriului lui Raabe – Duhamel aduce o simplificare a modului de calcul a limitei. 98 2 n 2 n 7.3 SERII ALTERNATE Definiţia 7.3.1. O serie u , u se numeşte serie alternată (cu termeni alternaţi) n 1 n n dacă, termenii săi satisfac relaţiile: un un 1 0 , ( )n (7.3.1) Observaţii: i) Din (7.3.1) rezultă că termenii seriei alternează ca semn; atunci termenul general un poate fi scris sub una din următoarele forme: un ( 1)n 1an, an 0 (dacă primul termen al seriei este pozitiv) respectiv: un ( 1)nan , an 0 (dacă primul termen al seriei este negativ) ii) Deci orice serie alternată poate avea una din următoarele două forme: ( n 1 ( n 1 1)n 1an a1 a2 a3 a4 1)n an a1 a2 a3 a4 ( 1)n 1an (7.3.2) ( 1)n an Exemple: a) ( 1)n 1 1 1 n 1 1 1 1 22 32 42 ( 1)n 1 1 b) ( 1)n n 1 2 3 4 ( 1)n n n 1 În studiul convergenţei seriilor alternate nu pot fi folosite criteriile de la serii cu termeni pozitivi. Vom prezenta aici un singur criteriu specific seriilor cu termeni alternanţi. 99 n 3 n Criteriul lui Leibniz: Fie seria alternată ( n 1 1)nan cu an 0. Dacă şirul (an)n satisface condiţiile: i) an ii)(an)n 0 estemonotondescrescător (an an 1, ( )n ) atunci seria ( n 1 1)n an converge. Observaţii: i) Criteriul lui Leibniz oferă numai condiţii suficiente de convergenţă nu şi necesare; cu alte cuvinte dacă condiţiile nu sunt satisfăcute nu rezultă că seria diverge, ci doar faptul că criteriul nu poate fi aplicat. ii) Totuşi, dacă lim a n 0, atunci ( n 1 1)n an (D), dar datorită aplicării criteriului general de divergenţă lim u lim( n n n 1)n an 0 . iii) În cazul în care an 0, dar şirul (an)n nu este monoton descrescător criteriul lui Leibniz nu este satisfăcut, şi nu putem afirma nimic despre natura seriei ( n 1 1)nan . Exemple: a) ( 1)n 1 n 1 1 n 1 1 1 1 2 3 4 ( 1)n 1 1 (seria armonică alternată) n Avem şirul an 1 n 0 şi este monoton descrescător, deci conform criteriului lui Leibniz seria converge. b) ( n 1 1)n 3nn 1 ; cu an n 3n 1 1 deci prima condiţie a criteriului lui Leibniz nu este satisfăcută, acesta neputând fi aplicat. Totuşi, deoarece: n lim u lim( n n n 1)n 3n+1 nu există, rezultă conform C.G.D. că seria diverge. c) ( n 1 1)n 1 sin2 n , cu an sin2 n n 0 100 n Deoarece 0 sin2 n 1 0 an 1 . Trecând la limită (teorema cleştelui) obţinem n lim a 0, deci prima condiţie din criteriul lui Leibniz este satisfăcută. Se poate arăta relativ n n uşor că şirul (an)n nu este monoton, deci condiţia a doua nu este satisfăcută şi criteriul lui Leibniz nu poate fi aplicat. Pentru a putea totuşi preciza natura acestei serii, se pot aplica criterii de la serii cu termeni oarecare (neprezentate aici) sau folosi noţiunile din următorul subcapitol. 7.4 SERII ABSOLUT CONVERGENTE. SERII SEMI – CONVERGENTE. Evident, în cazul seriilor cu termeni oarecare a , n 1 an , nu se pot aplica în studiul naturii acestora, criteriile formulate pentru seriile cu termeni pozitivi sau pentru seriile cu termeni altenanţi, enunţate în secţiunile precedente 7.2 respectiv 7.3. Totuşi, în secţiunea 7.1 au fost enunţate două rezultate care pot fi aplicate tuturor seriilor numerice, şi anume: Criteriul general de convergenţă a lui Cauchy, respectiv Criteriul general de divergenţă. Aplicarea lor, în special pentru stabilirea convergenţei seriilor, este aşa cum am mai spus, de cele mai multe ori extrem de dificilă. De aceea, matematicienii au căutat şi au stabilit, noi criterii, mai simple de aplicat în practică, cum ar fi: Criteriul lui Abel, Criteriul lui Dirichlet, etc. Chiar şi aşa, aplicarea lor necesită abilităţi matematice mai deosebite, precum şi o oarecare practică în lucrul cu serii. Deoarece lucrarea de faţă se adresează în mod special studenţilor economişti, s-a considerat utilă renunţarea prezentării acestor noi criterii; vom arăta totuşi că, în anumite situaţii, studiul convergenţei / divergenţei seriilor cu termeni oarecare se reduce la studiul convergenţei / divergenţei seriilor cu termeni pozitivi. 101 n n n n n Definiţia 7.4.1. Fie seria cu termeni oarecare a , n 1 an . Spunem că: a) Seria a , este absolut convergentă (A.C.), dacă seria moduleleor a este n 1 n 1 n n convergentă (ca serie cu termeni pozitivi). b) Seria a , este semi - convergentă (S.C.), dacă a (C) dar , a (D). n 1 n 1 n 1 n n n Teorema 7.4.1. Dacă seria a , n 1 an este absolut convergentă, atunci şi seria n 1an este convergentă. Demonstraţie. Deoarece, seria a (A.C.) n 1 a (C). Aplicând Teorema 7.1.3. n 1 seriei modulelor, obţinem: ( ) 0, ( )n a.î. an 1 an k an 1 an k , ( )n n şi ( )k (*) Din proprietăţile funcţiei modul, rezultă că: an 1 an k an 1 an k (**) Din (*) +(**) rezultă că a n 1 Cauchy. (C), conform criteriului general de convergenţă a lui q.e.d. Observaţie: Reciproca teoremei de mai sus nu este adevărată; cu alte cuvinte, convergenţa seriei simple nu implică convergenţa seriei modulelor a (C) a (C) . n 1 n 1 n n 102 n n n n n n n n n n Exemplu: Fie seria a n 1 ( 1)n 1 . Evident este o serie alternată şi converge deoarece n 1 satisface criteriul lui Leibniz. Totuşi, seria modulelor, a n 1 simplă. 1 , diverge ca serie armonică n 1 1. Întrucât seria valorilor absolute, a , este o serie cu termeni pozitivi, pentru studiul n 1 n naturii acesteia se vor folosi criteriile enunţate în secţiunea 4.3 2. În cazul în care seria modulelor a n 1 (D), am văzut că despre seria simplă a , nu se n 1 poate afirma nimic. Totuşi, se demonstrează că în cazul în care divergenţa seriei modulelor a a fost stabilită cu Criteriul raportului sau cu Criteriul rădăcinii, atunci n 1 n şi seria (simplă) a diverge. n 1 n 3. Pentru seriile cu termeni pozitivi a , n 1 an 0 , noţiunile de convergenţă şi convergenţă absolută coincid (deoarece a a ). n 1 n 1 n n 4. În anumite situaţii, este utilă şi următoarea formulare echivalentă a Teoremei 4.7. “Dacă seria a , n 1 an diverge, atunci şi seria a diverge”. n 1 În concluzie, în studiul naturii seriilor cu termeni oarecare, a , n 1 aplica următorul: Mod de lucru: an , se poate i) Se ataşează seriei cu termeni oarecare, seria modulelor a (care este serie cu n 1 n termeni pozitivi); 103 n n n ii) Se studiază natura seriei a folosind criteriile enunţate pentru serii cu termeni n 1 n pozitivi; iii) Dacă: a) a (C) a (C) n 1 n 1 n n b) a (D) (divergenţa fiind stabilită cu Criteriul raportului sau Criteriul n 1 n rădăcinii a (D) n 1 n c) n 1 an (D) (divergenţa fiind stabilită cu alte criterii) despre natura seriei a nu se poate spune nimic. În acest caz trebuie folosite alte criterii (Abel, n 1 n Dirichlet etc.). Exemple: 1) ( 1)n 2n sinn , cu 0,2 parametru real. n 1 n Evident seria este o serie cu termeni oarecare (datorită lui sinn ) şi nu cu termeni alternanţi. Ataşăm seria a n 1 2n sinn n 1 n şi aplicăm criteriul raportului. Atunci, obţinem: lim an 1 lim 2n 1 sinn 1 n n n a n 1 n 2n sinn lim 2 n n n 1 sin şi deci: 2 sin Atunci, pentru: i) 2 sin 1, seria a (C) n 1 a (C) (conf.T.7.1) n 1 104 n n n n n n n n , n ii) 2 sin 1, seria a (D) n 1 a (D) (conf.Obs.3) n 1 iii) 2 sin 1, nu putem preciza natura seriei a n 1 şi deci a seriei iniţiale a . n 1 Rezovând inecuaţiile trigonometrice de mai sus, avem: i) pentru 0, 6 2 7 11 , 3 6 6 ,2 ( 1)n 2n sinn (A.C.) n 1 ii) pentru 2 7 11 , , 6 3 6 6 ( 1)n 2n sinn (D) n 1 iii) pentru 2 7 11 6 3 6 6 , , nu putem preciza natura seriilor n 1 an şi n 1an . Dar, în cazul iii) cele două serii (simplă şi a modulelor) au expresiile: a) an ( 1)n n 1 n 1 1 n (C), ca serie alternantă care satisface criteriul lui Leibniz (şirul 1 n 0, monoton descrescător). b) a n 1 1 1 n 1 n n 1 n1/2 (D), (ca serie armonică 1 cu n 1 1 2 1) Deci, în acest caz, seria ( 1)n 2n sinn este (S.C.) n 1 n 2) n 1 n 1 a n n 1 n 1 2a cu a \ 1 2 parametrul real. Ataşăm seria modulelor: a n 1 n 1 n 1 n 1 n 1 n 2a şi aplicăm criteriul rădăcinii: a lim n a lim n n n n Atunci, pentru: n 1 n 1 a n n 1 2a lim nn 1 1 1 n a 1 a 2a 1 2a i) ii) 1 a 1 2a 1 a 1 2a 1 a (C),deci şi seria a (C) (conf. T.4.7) n 1 n 1 n n 1 a (D), deci şi seria a (D) (conf. Obs.3) n 1 n 1 n n 105 n iii) 1 a 1 2a 1 nu putem preciza natura seriilor a şi a . n 1 n 1 n n Atunci, rezolvând inecuaţiile de mai sus, concluzionăm că: i) pentru ,0 2, 3 n 1 n 1 n 1 n 1 n 2a este (A.C.); a ii) pentru 0, 2 \ 3 1 n 1 n 1 2 n 1 n 1 n 2a este (D); a iii) pentru 0, 2 nu putem preciza natura seriei. 3 Dar, pentru: a) 0 a n 1 n 1 n , care (D) conform C.G.D. n 1 n Într-adevăr lim a lim n n n n 1 n n lim 1 n 1 n n e 0 b) 2 3 a ( 1)n n 1 n 1 n n 1 n n , care (D) conform C.G.D., deoarece an 0. 106 US8 Serii de puteri Seriile de puteri sunt un caz particular de sume infinite. Ele sunt utilizate pentru aproximarea funcțiilor transcendente. Capitolul de față abordează modul în care poate fi studiată o serie de puteri: determinarea razei de convergență, studiul convergenței în extremitățile intervalului de bază, determinarea sumei, proprietățile sumei. Studentul trebuie să rețină metoda de aflare a razei de convergență și să realizeze că în capetele intervalului de bază obținem o serie numerică de tipul celor studiate în US7. Pentru partea teoretică două ore sunt suficiente. Pentru aplicații sunt necesare măcar patru ore: trei pentru mulțimea de convergență și una pentru studiul in capetele intervalului. Utile din lista bibliografică se dovedesc a fi lucrarile de la pozițiile [3], [6] și [9]. * * * Aşa cum s-a observat din secţiunile precedente, noţiunea de “serie numerică” apare ca o generalizare (extindere) firească a noţiunii de “şir numeric”. În mod natural, aceste două noţiuni au căpătat în Analiza matematică noi extinderi, din ce în ce mai complexe, lucru de altfel deseori întâlnit nu numai în matematică ci şi în multe alte domenii ale ştiinţei. Astfel, pornind de la noţiunea de şir numeric: (an)n :a0,a1,a2, ,an (*) se introduce noţiunea de “şir de funcţii”: fn(x) n :f0(x),f1(x), ,fn(x), (**) unde funcţiile fn (x) sunt definite pe acelaşi domeniu D . Exemple: a) fn(x) n , cu termenul general: fn : fn(x) nx2n sinnx adică, şirul de funcţii are forma: xsinx, 2x4 sin2x, 3x6 sin3x, ,nx2n sinnx, 107 n n b) fn(x) n , unde termenul general este: fn :(0, fn(x) ) (n2 1)ln(nx) deci, şirul de funcţii are forma: 2lnx, 5ln2x, 10ln3x, ,(n2 1)ln(nx) Ţinând seama că noţiunea de “funcţie “ este mult mai complexă decât cea de “număr”, este evident că noţiuni precum convergenţă sau divergenţă au căpătat pentru şirurile de funcţii valenţe noi. Astfel, un şir de funcţii fn(x) n poate fi: simplu (punctual) convergent, uniform convergent, slab convergent, tare convergent etc. Este natural ca în caz de convergentă, şirul de funcţii să aibă limita o funcţie f(x) şi nu un număr. De asemenea, ţinând cont că termenii şirului sunt funcţii, care pot avea diverse proprietăţi (funcţie continue, derivabile, integrabile etc.), este natural să cercetăm în ce condiţii acestea se “transferă” la funcţia limită f(x) . În mod cu totul analog noţiunea de “serie numerică”: a n 0 a0 a1 a2 an (***) capătă o primă generalizare şi anume cea de “serie de funcţii”: f (x) n 0 f0(x) f1(x) f2(x) fn(x) (****) unde termenii seriei sunt funcţii definite pe acelaşi domeniu D . Exemple: a) (n2 1)arctg(xn 1) arctg1 2arctg(x 1) 5arctg(x2 1) unde termenul n 0 general al seriei este: fn : fn(x) (n2 1)arctg(xn 1) b) n! 1! 2! 3! n 1 1 xn 1 x 1 x2 1 x3 108 n n n 0 f : 0,1 n unde termenul general are forma: fn (x) 1 nxn ! Analog, în caz de convergenţă (simplă, uniformă, etc.) seria va converge la o “funcţie sumă” s(x). Din nou este interesant şi important de studiat, în ce condiţii proprietăţile termenilor seriei de funcţii (integrabilitate, derivabilitate, continuitate etc.) se transferă asupra funcţiei sumă s(x). Un caz particular de serii de funcţii, dar foarte important şi des intâlnit în aplicaţiile economice, este cel al seriilor de puteri, unde termenii seriei sunt funcţii putere. Astfel, avem: Definiţia 8.1: Numim serie de puteri centrată în x0 , o serie de funcţii de forma: n 0an(x-x0)n a0 a1(x-x0) a2(x-x0)2 an(x-x0)n (8.1) unde (an)n este un şir numeric (Numerele an se numesc coeficienţii seriei de puteri). Observaţii: i) un caz particular, dar foarte des întâlnit în aplicaţii, este cel al seriilor de puteri centrate în 0 (pe scurt, serii de puteri) deci când x0 0. Vom obţine atunci seria: n 0anxn a0 a1x a2x2 anxn (8.2) ii) Deoarece prin substituţia: y x x0 (8.3) o serie de puteri centrată în x0 (şi variabila x) se reduce întotdeauna la o serie de puteri centrată în 0 (şi variabila y), adică avem: a (x n 0 x0)n a yn (8.4) n 0 ne vom ocupa numai în continuare numai de seriile de puteri centrate în 0, deci de forma (8.2). Definiţia 8.2: Spunem că, seria de puteri a xn este: n 0 n a) convergentă (punctual) în x x0 convergentă; , dacă seria numerică ataşată, a xn este n 0 109 n n b) convergentă (punctual) pe domeniul D0 , dacă converge în orice punct x D0 (D0 se numeşte mulţime de convergenţă a seriei). Teorema 8.1. Mulţimea de convergenţă a unei serii de puteri, a xn , este n 0 n întotdeauna nevidă. Demonstraţie. Fie seria de puteri n 0anxn cu (an)n un şir de numere reale oarecare şi D0 mulţimea sa de convergenţă. Pentru x deci D0 0 seria a xn n 0 . a0 , deci converge şi are suma S a0 . Rezultă că 0 D0 , q.e.d. Teorema 8.2. (lui Abel) Fie seria de puteri, a xn , n 0 convergenţă), astfel încât: an . Atunci, ( !)r 0, (“r” se numeşte rază de i) pentru x ii) pentru x ( r, r) a xn converge (C) n 0 (8.5) n ( , r) (r, ) a xn diverge (D) n 0 n Observaţii: i) Teorema lui Abel nu precizează natura seriei în capetele intervalului de convergenţă: (-r, r). În exemplele concrete, se studiază separat natura seriei în cele două capete, adică seriile numerice: a rn şi a ( r)n . n 0 n 0 n n 110 n n ii) Dacă raza de convergenţă r 0, atunci evident seria de puteri a xn converge n 0 n numai în x 0 şi diverge în rest (pentru x ). iii) Dacă raza de convergenţă r , atunci mulţimea de convergenţă a seriei D0 , deci seria nu diverge pentru nici o valoare reală. iv) Aşadar, folosind teorema lui Abel, dacă aflăm raza de convergenţă r, ştim mulţimile pe care seria a xn converge sau diverge (evident cu excepţia valorilor x r ). r 0 n Următoarele două teoreme ne oferă două formule pentru aflarea razei de convergenţă r. Teoremă 8.3 (Teorema raportului). Fie seria de puteri, a xn cu a n 0 n convergenţă r , este dată de relaţiile: şi lim an 1 0, . Atunci raza de n n a 1 ; (0, ) r ; 0 (8.6) 0 ; Teorema 8.4 (Teorema rădăcinii) Fie seria de puteri, a xn n 0 cu an , şi lim n a 0, . Atunci, raza de n n convergenţă r, este dată de relaţiile: 111 0 n n n n 1 ; (0, ) r ; 0 (8.7) 0 ; Observaţie: Din rezultatele enunţate anterior, se observă că în determinarea mulţimii (intervalului) de convergenţă respectiv a mulţimii (reuniune de intervale) de divergenţă pentru seria de puteri a xn , trebuie aplicat următorul: n 0 n Algoritm de lucru: 1) Se calculează valoarea “ ” aplicând una din formulele: lim an 1 sau lim n a a n n n n 2) Să determină raza de convergenţă “r” cu relaţiile (8.6) sau (8.7) (care de fapt coincid); 3) Se scriu intervalele de convergenţă / divergenţă ale seriei a xn ; n 0 n 4) Se studiază natura seriei în cele două capete x r . Observaţie: În cazul în care seria este centrată în x a (x n 0 x0)n se reduce la seria centrată în 0 a yn folosind substituţia (8.3) şi apoi se aplică algoritm de lucru, obţinându-se n 0 n mulţimea de (C)/(D) pentru variabila y. Folosind aceeaşi substituţie se scriu mulţimile de (C)/(D) în variabila x. Exemple: a) ( n 1 1)n xn ; Avem: an ( 1)n 1 , deci an 1 . Folosind teorema rădăcinii, obţinem: 112 n n n n n n n a lim n a lim n 1 n n lim n1 1 (am folosit aici limita fundamentală: lim n n 1) n deci conform (8.7), avem: r 1 1 deci r 1. Atunci, conform Teoremei 8.2. (Abel) 1 obţinem: i) pentru x (1, 1) ( 1)n xn (C) n 1 n ii) pentru x ( , 1) (1, ) ( 1)n xn (D) n 1 n iii) pentru x 1 ( 1)n xn nu-i putem preciza natura. n 1 n Dar, pentru: x 1 obţinem: a xn n 1 ( 1)n 1 (C) n 1 (conform Criteriului lui Leibniz de la serii alternante) iar, pentru: x 1, obţinem: a xn n 1 1 n 1 n (D)(ca serie armonică n 1 1 n cu =1) În concluzie, ( 1)n xn n 1 n (C), dacă x (D), dacă x ( 1,1] ( , 1] (1, ) 2) ( n 0 n 1)n 1 n! ; avem an ( 1)n 1 n! , deci an x 1 1 n! Folosind teorema raportului, obţinem: lim an 1 n n lim (nn 1)! lim n1 1 0 ! n n şi deci conform (4.6.6) avem r . Deci, seria ( n 0 n 1)n 1 n! (C) pentru ( ) x . x 113 n n n n a y n n n n n n 3) n 0 n 21(x 1)n , este o serie de puteri centrată în x0 1. Notăm y Atunci: x 1 şi obţinem n 0 n 21 yn cu an n 1 n 2n an lim an 1 n n lim (n n1)2n 1 2 n n2n 1 n2 2n 1 n 1 2 n n2 2n 1 2 lim Atunci raza de convergenţă r 1 2, şi avem: i) pentru y ( 2,2) n 1 n 1 n n2n (C) ii) pentru y ( , 2) (2, ) n 1 n 1 n2n (D) iii) pentru y 2 ?(nu putem preciza natura seriei) Revenind la variabila x(x y 1), avem: i) pentru x ( 3,1) n 1 n 21 (x+1)n (C) ii) pentru x ( , 3) (1, ) n 1 n 21 (x+1)n (D) iii) pentru x 3, x 1 ?(nu putem preciza natura seriei) Dar, pentru: x 3, obţinem a (x 1)n ( n 1 n 1 1) nn 1 (D)(conform C.G.D, deoarece termenul general ( 1)n nn 1 0) iar pentru: x = 1, obţinem a (x 1)n n 1 n 1 n 1 n (D) (conform C.G.D, deoarece termenul general nn 1 1 0). 114 US9 Funcții de mai multe variabile Funcțiile de mai multe variabile constituie generalizarea firească a funcțiilor de o variabilă reală. Studiul lor se justifică prin prezența acestora în majoritatea modelelor economice. În prezenta unitate de studiu se extind treptat noțiuni familiare studenților: domeniu de definiție, grafic, limită, derivată, diferențială, punct de extrem. Am presupus cunoscute doar noțiunile fundamentale de analiză matematică și am prezentat aproape totul în material. Credem că astfel ajutăm studenții proveniți din licee care au în programă mai puțină matematică să își însușească și să aplice aceste metode. Noțiunea centrală este cea de punct de extrem cu care capitolul se încheie. Studenții trebuie în principiu să învețe să privească funcțiile de mai multe variabile ca o modalitate mai riguroasă de a descrie realitatea. Ei trebuie să învețe să calculeze derivatele parțiale (de diferite ordine) ale unei astfel de funcții, să determine diferențialele (de ordinul I și II) și să le aprecieze rolul în studiul funcțiilor, să determine și să caracterizeze punctele de extrem. Este un capitol amplu, cu multe noțiuni noi și cu nenumărate aplicații. Pentru acest capitol zece ore pot fi suficiente: două ore pentru partea teoretică, două ore pentru deprinderea calculului derivatelor parțiale de ordinul I și a diferențialei, două ore pentru derivatele parțiale de ordin superior și a diferențialei a 2-a, patru ore pentru aplicațiile la punctele de extrem. Orice manual universitar de analiză matematică are aceste noțiuni incluse. Pentru doritori trimiterile bibliografice [1], [3] și [6] pot fi de asemenea folositoare. * * * Fenomenele vieţii reale sunt descrise cel mai adesea de funcţii care depind nu de o singură variabilă ci de mai multe variabile independente. Acest lucru este valabil şi în viaţa economică. În acest capitol se introduc elementele de bază ale teoriei acestui tip de funcţii. Pentru a depăşi dificultăţile extinderii conceptelor fundamentale (continuitate, limită, derivabilitate etc.) de la funcţii de o variabilă la o funcţie de mai multe variabile am abordat numai cazul funcţiilor de două variabile, dezvoltarea ulterioară făcându-se în aceeaşi manieră. 115 9.1. ŞIRURI DE PUNCTE ÎN n Să considerăm spaţiul vectorial n şi să considerăm un vector oarecare x (x1,x2, ,xn)t al său. Pentru a simplifica scrierea vom renunţa la simbolul de transpunere şi în cele ce urmează vom scrie pur şi simplu x (x1,x2, ,xn). Acestui element al spaţiului n îi asociem un punct P care are coordonatele x1,x2, ,xn , tot aşa cum unui vector x (x1,x2) din 2 îi asociem un punct P cu coordonatele x1,x2 într-un sistem de axe ortogonale. Asociem unei perechi de astfel de puncte P, P0 P(x , x , , x, )P (x1 ,x2, ,x0) un 1 2 0 0 0 n n număr nenegativ numit distanţă şi definit prin: d(P,P ) (x1 x0)2 (x2 x0)2 (xn x0)2 . (9.1.1) 0 1 2 n Evident, distanţa este nulă dacă şi numai dacă P şi P0 coincid. Convenim că însuşi vectorul x să se numească punct din n . În acest mod, distanţa d se asociază unei perechi de vectori din n . Definiţia 9.1.1. Se numeşte şir de puncte din n funcţia care asociază oricărui număr natural k un punct unic xk din n . k xk (xk,xk, ,xk) 1 2 n Şirurile (x1 )k , (x2)k , ,(xn)k se numesc şiruri de coordonate. k k k Definiţia 9.1.2. Şirul (xk)k converge la punctul x0 dacă: lim d(xk,x0) 0 (9.1.2) k În spiritul definiţiei cunoscute a limitei unui şir de numere reale aceasta înseamnă: ( ) 0, k( ) a.î.( )k k( ) d(xk,x0) (9.1.3) Folosind (9.1.1) deducem că: (xk x0)2 (xk x0)2 (xk x0)2 1 1 2 2 n n 116 1 , 2 3 pentru k k( ) adică xk x0 ( ) k k( ), ( ) i 1,n (9.1.4) i i sântem conduşi astfel la: Teorema 9.1.1. Şirul xk este convergent la x0 dacă şirurile de coordonate xi sunt k convergente la x0 pentru orice i 1,n. i Exemplu. Considerăm în 2 şirul xk 2k 1 k 2 k2 3k , . Şirurile de coordonate sunt xk 2k 1 k2 xk k3k2 şi pentru k au limitele 0 şi respectiv 1 . Aceasta înseamnă că xk are limita x0 0,1 . 3 9.2. FUNCŢII DE DOUĂ VARIABILE În cele ce urmează vom considera spaţiul 2 ale cărui puncte x au doar două coordonate x (x1,x2). Pentru ca scrierea să fie mai simplă vom înlocui perechea de coordonate (x1,x2 ) prin perechea (x,y). Fie D o mulţime de puncte (x,y) din 2 şi I o mulţime din . Definiţia 9.2.1. Se numeşte funcţie de două variabile tripleta formată din D,I şi o lege de corespondenţă care asociază oricărui punct (x,y) D, un singur element z I. Exprimăm această corespondenţă prin z f(x,y). Observaţie. Când D este din n iar I în mod similar se defineşte o funcţie de n variabile. Dacă luăm un sistem cartezian de coordonate Oxyz şi reprezentăm mulţimea punctelor P(x,y,z) unde z f(x,y) obţinem imaginea geometrică a funcţiei sau graficul său. De data aceasta graficul e o porţiune dintr-o suprafaţă. 117 Fig. 5.1. Exemplu: Să se reprezinte graficul funcţiei z x2 y2 . Soluţie. Ridicăm la pătrat egalitatea şi obţinem z2 x2 y2 . Dacă facem aici x = 0 (intersectăm graficul cu planul yOz obţinem z y adică două drepte în acest plan. În mod similar intersectăm graficul cu planul xOz şi obţinem dreptele x z. Graficul este un con infinit cu vârful în originea sistemului de coordonate. Fig. 9.2 118 y x x x 9.3. LIMITĂ ŞI CONTINUITATE Să considerăm o funcţie f :D 2 şi (x0,y0) un punct de acumulare al lui D. Aceasta înseamnă că există cel puţin un şir de puncte (xk,yk)k din D care converge la (x0,y0). Definiţia 9.3.1. Funcţia f(x,y) are în (x0,y0) limita , dacă oricare ar fi şirul (xk,yk)k D, (xk,yk)k (x0,y0), (xk,yk) (x0,y0) rezultă că f(xk,yk) . Vom scrie: lim f(x,y) (x,y) (x0,y0) sau lim f(x,y) . x x y y0 0 Această limită se numeşte limită globală deoarece x şi y tind simultan şi independent către x0 respectiv y0. Dacă x0 şi y tind succesiv la x0 respectiv y0 adică dacă 1 lim ( lim f(x,y)), x x0 y y0 2 lim ( lim f(x,y)),, se obţin limitele iterate sau parţiale. Are loc: y y0 x x0 Teorema 9.3.1. Dacă limita globală există atunci şi limitele parţiale există şi are loc egalitatea: 1 2 (9.3.1) Reciproca nu este adevărată. Demonstraţie. Pentru a justifica rezultatul este suficient să considerăm mai întâi şiruri generale (xk,yk)k (x0,y0) şi apoi să le particularizăm în şiruri de formă particulară (xk,y0)k (x0,y0) şi (x0,yk)k (x0,y0). Ceea ce e adevărat pentru general se păstrează şi pentru particular. Faptul că reciproca e falsă se poate proba prin următorul exemplu. Să se calculeze limitele iterate ale funcţiei f(x,y) x x y în punctul (0,0). 1 lim lim x x 0 y 0 y y lim x lim1 1 x 0 x 0 2 lim lim x y 0 x 0 y y lim y 0 y y lim( 1) 1 y 0 Cele două limite nefiind egale nu există limita globală a funcţiei în origine. 119 k k , Observaţii: 10. Dacă 1 2 nu rezultă că limita globală există ci doar că dacă aceasta există este egală cu valoarea lor comună. 20. S-ar putea ca în punctul (x0,y0) funcţia f să nu fie definită adică să nu se poată calcula f (x0,y0) cum este în exemplul precedent. Definiţia 9.3.2. Funcţia f (x,y) este continuă în (x0,y0) dacă: lim f(x,y) x x y y0 0 f(x0,y0) Ca şi în cazul limitei globale acest tip de continuitate se numeşte continuitatea globală. Se poate şi aici pune problema continuităţii în raport cu fiecare argument când celălalt este fixat, adică continuitatea parţială. Exemplu. Să se studieze continuitatea globală şi cea parţială pentru funcţia: f(x,y) x2xyy2 ,(x,y) (0,0), 0 ,(x,y) (0,0) Calculăm limitele parţiale: 1 2 limf(x,0) lim0 0 x 0 x 0 limf(0,y) lim0 0 y 0 y 0 Se observă că 1 f(0,0) 2 0, deci funcţia este continuă în raport cu fiecare argument. Limita globală nu există pentru că luând şiruri (xk,yk)k (0,0) în care yk mxk, m 0 obţinem: f(xk,yk) xk mxk m x2 m2x2 1 m2 deci f(xk,yk) m 1 m2 care depinde de m. Aşadar continuitatea globală nu există. 120 y 9.4. DERIVATE PARŢIALE ŞI DIFERENŢIALE 9.4.1. DERIVATE Să considerăm funcţia f :D PARŢIALE DE ORDINUL ÎNTÂI şi (x0, y0) un punct din interiorul mulţimii D. Considerăm de asemenea un punct (x, y) din D diferit de (x0, y0). Definiţia 9.4.1.1. Se numeşte derivată parţială a funcţiei f în raport cu variabila x în punctul (x0, y0). lim f(x,y0) x x x0 f(x0,y0) x0 (9.4.1) dacă limita există. Dacă aceasta este şi finită, spunem că funcţia este derivabilă parţial în raport cu x în punctul considerat. Notăm valoarea limitei cu f (x0,y0) sau fx(x0,y0). În mod similar se introduce x derivata parţială şi derivabilitatea parţială în raport cu y în (x0, y0): f (x0,y0) = fy(x0,y0) lim f(x0,y) y y0 y f(x0,y0) y0 Observaţie. Derivata parţială nu aduce nici o noutate faţă de cazul funcţiilor de o singură variabilă. Practic, când se derivează parţial o funcţie de n variabile în raport cu una, toate celelelte variabile se consideră constante (parametri). Dacă f (x, y) este derivabilă parţial în raport cu x sau/şi y în toate punctele unei submulţimi D1 a a lui D, spunem că f este derivabilă parţial pe D1. În acest caz pentru orice punct (x0,y0) D1 stabilim corespondenţele (x0,y0) f (x0,y0), (x0,y0) f (x0,y0) care definesc două funcţii noi numite derivatele x y parţiale ale lui f în raport cu x respectiv y în D1. Ele se vor nota: f (x,y) sau fx(x,y) respectiv f (x,y) sau fy(x,y). x y 121 y x Observaţie. Calculul derivatelor parţiale f (x,y) , f (x,y)se face folosind regulile x y obişnuite de calcul pentru derivata sumei, produsului, câtului, funcţiilor compuse etc. Exemplu. Să se calculeze f (x,y) , f (x,y) dacă f(x,y) xy exy . x y Soluţii. Calculăm f (x,y) considerându-l pe y parametru constant. x f (x,y) =(xy)xexy (xy)(exy)x yexy xy(yexy) (y xy2)exy x Calculăm f (x,y) considerându-l pe x parametrul constant: y f (x,y) =(xy)yexy (xy)(exy)y xexy (xy)(xexy) (x x2y)exy y Observaţii : 10. În exemplul de mai sus, simetria în raport cu variabilele x şi y permitea calculul celei de-a doua derivate parţiale schimbând între ele variabilele x cu y în f (x,y). x 20. Toate consideraţiile anterioare se pot extinde de la funcţii de două variabile la funcţii de n variabile fără nici un fel de dificultate. 9.4.2. DERIVATE PARŢIALE DE ORDIN SUPERIOR Să considerăm din nou f :D 2 , derivabilă parţial în raport cu x şi y pe D. Cele două derivate parţiale fx (x,y)  g1(x,y) şi fy (x,y)  g2(x,y) pot fi la rândul lor derivabile pe   o submulţime D1 a lui D. Apare natural să considerăm derivatele lor parţiale drept derivate parţiale de ordinul al doilea ale funcţiei f: g1 (x,y) x f (x,y) 2f (x,y) fxx(x,y), x x 2 x g1 (x,y) y f (x,y) 2 y x (x,y) fxy(x,y), f 122 y y x y g2 (x,y) x f (x,y) 2 x y(x,y) fyx(x,y), f g2 (x,y)   f (x,y)  yf (x,y)  fy(x,y). 2     y y y y 2  Pentru funcţia de două variabile f(x,y) obţinem astfel patru derivate de ordinul al doilea. Dintre acestea se remarcă f (x,y) şi fx(x,y) numite derivate parţiale mixte de ordinul al xy y doilea. Natural se pune întrebarea dacă între cele două derivate există legături. Exemplul următor sugerează că răspunsul este afirmativ. Exemplu. Să se determine derivatele parţiale de ordinul al doilea ale funcţiei f :D 2 , f(x,y)  x2y2xy 1. Soluţie. Obţinem succesiv: f  2xy  2y ; f  x2  2x; x x 2f  xx  2y ; xy  xy  2x  2 2        f  f  f 2 x   yx  yx 2x  2 ; 2f  yy  0. 2 2   y      f  f  f   Se observă că cele două derivate parţiale mixte f ,fx sunt egale. Nu orice funcţie f xy y (x,y) are această proprietate. Fără demonstraţie enunţăm următoarea teoremă: Teorema 9.4.2(Criteriul lui Schwarz) Dacă în domeniul D1  D funcţia f (x,y) satisface condiţiile: a) are derivatele f , f , b) derivatele mixte de ordinul al doilea 2f 2f xy yx , sunt continue global, atunci derivatele parţiale mixte de ordinul al doilea sunt egale. 123 v Observaţie. Funcţia f (x,y) din exemplul anterior satisface condiţiile criteriului lui Schwarz. Fiecare derivată parţială de ordinul al doilea poate fi la rândul său derivată parţial în raport cu cele două variabile. Se obţin astfel derivate parţiale de ordinul al treilea şi procesul poate continua. 9.4.3. DERIVATELE PARŢIALE ALE FUNCŢIILOR COMPUSE Să luăm funcţia f :D 2 şi să presupunem că la rândul lor, argumentele x şi y ale lui f sunt funcţii de variabilele u,v,w: f(x,y)  f(x(u,v,w), y(u,v,w))  (u,v,w) (9.4.2) Ne interesează derivatele parţiale ale funcţiei  în raport cu cele trei variabile u,v,w. Un calcul nu foarte complicat care foloseşte definiţia derivatei parţiale conduce la formula de calcul a celor trei derivate parţiale:  f x f y u x u y u     , f x f v x v y f x f w x w y y , y . w Exemplu. Să se afle derivatele parţiale ale funcţiei f(x,y)  xy2  x2y, cunoscând că x = u + v, y = u - v. Soluţie. Notăm (u,v)  f(x(u,v), y(u,v))   f  x  f  y  (y2 2xy)1(2xy  x2)1 y2  4xy  y2 , u x u y u   f  x  f  y  (y2  2xy)1(2xy  x2)1 y2  x2 . v x u y v 124 x y Observaţie. Formula (9.4.3.1) se poate generaliza când f are n variabile iar fiecare este funcţie de m variabile. Evident aplicarea formulei este posibilă când toate funcţiile implicate sunt derivabile parţial. 9.4.4. DIFERENŢIALELE FUNCŢIILOR DE MAI MULTE VARIABILE Fie f :D 2 , (x0,y0) şi (x,y) două puncte din D. Să notăm h = x - x0, k = y - y0. Numim creşterea totală a funcţiei în punctul (x0,y0) cantitatea f(x0,y0) definită prin: f(x0,y0)  f(x,y)f(x0,y0) (9.4.3) Definiţia 9.4.2 Funcţia f este diferenţiabilă în (x0,y0) dacă există două constante A şi B independente de h şi k şi două funcţii (x0,y0,h,k),(x0,y0,h,k) astfel încât: (x0,y0,0,0)  0 , (x0,y0,0,0)  0, lim(x0,y0,h,k)  lim(x0,y0,h,k)  0 k0 k0 h0 h0 f(x0,y0)  AhBhh(x0,y0,h,k)k(x0,y0,h,k) (9.4.4) Se poate demonstra următoarea teoremă: Teorema 9.4.3 Dacă f (x,y) este diferenţiabilă în (x0,y0) atunci ea este derivabilă parţial în raport cu x şi y în acest punct şi f (x0,y0) A, f (x0,y0) B. x y Definiţia 9.4.3 Expresia Ah +Bk se numeşte diferenţiala funcţiei f în punctul (x0,y0) şi se notează: df (x0,y0) = Ah + Bk. (9.4.5) Dacă funcţia f (x,y) este diferenţiabilă în (x0,y0) atunci diferenţiala sa df(x0,y0) e dată de df(x0,y0)  f (x0,y0)h  f (x0,y0)k Cum dx = h, dy = k obţinem diferenţiala de ordinul întâi finală: (9.4.6) a funcţiei f (x,y) în forma 125 2 x 2 y 2 x 2 y df(x0,y0)  f (x0,y0)dx f (x0,y0)dy x y Exemplu. Să se determine diferenţiala de ordinul întâi a funcţiei f(x,y)  xy2  2x3y în punctul (1,1). Soluţie: fx(x,y)  y2 6x2y  f(x,y)  2xy 2x3 y f(1,1)  7 x f(1,1)  4 y Obţinem: df(1,1) 7dx 4dy. Observaţie. Diferenţiala de ordinul întâi ale funcţiei f (x,y) serveşte la aproximarea creşterii totale a funcţiei f, f(x0,y0) conform formulei (9.4.3). Dacă f (x,y) este diferenţiabilă în toate punctele (x0,y0) unui domeniu D1 D, putem defini o corespondenţă: (x0,y0) df(x0,y0) ( ) (x0,y0) D1. Care se numeşte diferenţiala lui f şi se scrie df(x,y). Ea poate fi la rândul său diferenţiabilă. Diferenţiala sa se numeşte diferenţiala de ordinul al doilea a lui f şi se notează d2f(x,y). Un calcul elementar conduce la formula d2f(x,y). d2f(x0,y0) 2f (x0,y0)h2 2 2 x y(x0,y0)hk f 2f (x0,y0) k2 Exemplu. Să se determine d2f(1,1) pentru funcţia: f(x,y) xy 2x2y3 Soluţie. f x y 4xy3 f y x 6x2y2 2f x2 4y3 , 2f x y 1 12xy2 , 2f y2 12x2y 2f (1,1) 4 , 2 x y(1,1) 13 , f 2f (1,1) 12 şi în final 126 d2f(1,1) 4h2 26hk 12k2 . Observaţii: 10. Procesul poate continua obţinând diferenţiale de ordin superior lui 2. 20. Diferenţiala de ordinul al doilea este după cum se poate vedea din exemplul de mai sus o formă pătratică de variabile h,k. Când variabilele sunt x1,x2, ,xn această diferenţială devine o formă pătratică de n variabile. Ea joacă rolul derivatei a doua la funcţii de o singură variabilă adică serveşte studiul convexităţii unei funcţii de n variabile. 9.5. EXTREMELE FUNCŢIILOR DE MAI MULTE VARIABILE 9.5.1. FORMULA LUI TAYLOR PENTRU FUNCŢII DE MAI MULTE VARIABILE Se cunoaşte formula lui Taylor pentru funcţii de o variabilă care are derivate până la ordinul k+1 într-o vecinătate a punctului x0. f(x) f(x0) x 1x0 f (x0) (x 2x0)2 f (x0) k 1 ! ! f (x ) f ( ) (x x0)k (k) (x x0) (k 1) k! 0 (k 1)! fiind un punct intermediar lui x0 şi x. Această formulă permite aproximarea valorii f în x cu ajutorul valorilor acesteia şi a primelor sale k derivate în punctul x0, dacă ultimul termen (restul de ordin k) este neglijat. Căutăm o exitindere a formulei la cazul funcţiilor de două variabile după care, natural, extinderea se poate face la n variabile. Fie aşadar, f :D 2 şi (x0,y0) un punct interior din D. Presupunem că (x,y) este un punct oarecare din D şi notăm h x x0, k y y0 . 127 Teorema 9.5.1. Dacă f (x,y) are toate derivatele parţiale până la ordinul m în (x0,y0) iar derivatele parţiale de ordin m+1 există într-o vecinătate V a lui (x0,y0) atunci oricare ar fi punctul (x,y) V, are loc formula: f(x,y) f(x0,y0) 1! x (x0,y0)h y(x0,y0)k (2) 1 f f 1 2! x (x0,y0)h y (x0,y0)k f (9.5.1) (m 1) 1 m! x (x0,y0)h y (x0,y0)k f Rm unde : (i) x (x0,y0)h x (x0,y0)k f reprezintă puterea formală “i” a expresiei diferenţiale din paranteză aplicată funcţiei f iar exprimat prin: Rm numit restul de ordin m în formula lui Taylor este (m 1) Rm x (x0 h,y0 k)h y(x0 h,y0 k)k f; (0,1) (9.5.2) 9.5.2. TEOREMA LUI FERMAT PENTRU FUNCŢII DE MAI MULTE VARIABILE În teoria diferenţială a funcţiilor de o variabilă, teorema lui Fermat afirmă că într-un punct de extrem local din interiorul domeniului de derivabilitate al funcţiei, Df , derivata f se anulează. Încercăm aici o extindere la funcţii de două variabile care se poate apoi transpune la funcţii de n variabile. Să considerăm f :D 2 şi (x0,y0) D. Definiţia 9.5.1 Punctul (x0,y0) se numeşte maxim (minim) local pentru f(x,y) dacă există o vecinătate V a acestuia astfel încât ( ) (x,y) V are loc: f(x,y) f(x0,y0) (f(x,y) f(x0,y0)) (9.5.3) 128 Teorema lui Fermat poate constitui o metodă prin care se determină punctele de extrem sau posibilele puncte de extrem. Teorema 9.5.2 Dacă f(x,y) are derivate parţiale finite în raport cu x şi y în punctul de extrem local (x0,y0) din interiorul lui D, atunci: f (x0,y0) f (x0,y0) 0 (9.5.4) x y Demonstraţie. Considerăm (x) f(x,y0) şi (y) f(x0,y). Cum (x0,y0) este în interiorul lui D, (vezi fig. 5.3) x0 se află în interiorul lui (a,b) iar y0 în interiorul lui (c,d). Se poate aplica teorema lui Fermat pentru (x) cu x0 punct de extrem local şi (y) cu y0 punct de extrem local. Obţinem: dx (x0) f (x0,y0) 0, dy (y0) f (x0,y0) 0 d x d y ceea ce demonstrează teorema. Fig. 5.3 Soluţiile sistemului (9.5.4) se numesc puncte critice sau puncte staţionare. Reciproca nefiind adevărată, punctele critice rămân posibile puncte de extrem. Condiţiile (9.5.4) pentru punctele (x0,y0) sunt doar condiţii necesare de extrem. 129 Exemplu. Să se determine posibilele puncte de extrem ale funcţiei: f(x,y) x3 y2 3xy 15 Soluţie: Calculăm f , f . Avem: x y f x 3x2 3y, f y 2y 3x Rezolvăm sistemul f x 0, f y 0 3x2 3y 0 2y 3x 0 Obţinem soluţiile (0,0) şi 9 27 4 8 , care sunt posibile puncte de extrem. 9.5.3. CONDIŢII SUFICIENTE DE EXTREM (DETERMINAREA EXTREMELOR) Am văzut în cele de mai sus că un punct de extrem este punct staţionar dar din teorema lui Fermat nu putem obţine şi reciproca. Căutăm aşadar condiţii suficiente de extrem prin care să putem aprecia dacă un punct staţionar este sau nu este punct de extrem. Să remarcăm mai întâi că formula lui Taylor (9.5.1.2) se poate scrie şi sub forma: f(x,y) f(x0,y0) 1!df(x0,y0) 2!d2f m!dmf(x0,y0) Rm (9.5.5) Dacă m = 2 atunci 1 1 1 f(x,y) f(x0,y0) 1!df 2!d2f R2 (9.5.6) 1 1 Să presupunem acum că (x0,y0) este un extrem local. În virtutea condiţiilor (9.5.2.2) rezultă că : df(x0,y0) f (x0,y0)h f (x0,y0)k 0. x y 130 x 2 În plus dacă x x0 h şi y y0 k sunt suficient de mici, R2 poate fi neglijat (acolo h şi k apar la puterea a treia). Rezultă că pentru (x,y) într-o vecinătate suficient de mică a lui (x0,y0): sgn f(x,y) f(x0,y0) sgnd2f(x0,y0) obţinem, astfel: Teorema 9.5.3. Dacă forma pătratică d2f(x0,y0) este pozitiv definită în punctul critic (x0,y0) atunci acest punct este minim local. Dacă în acelaşi punct d2f(x0,y0) este negativ definită, punctul este maxim local. Condiţiile în care o formă pătratică este pozitiv (negativ) definită au fost prezentate anterior. Pentru prezenta formă pătratică, matricea este: f2 (x0,y0) H x fy (x0,y0) 2 x y(x0,y0) f 2 f y2 (x0,y0) (9.5.7) care se numeşte matrice hessiana şi se poate generaliza la funcţii de n variabile. Suntem conduşi astfel la următoarea teoremă: Teorema 9.5.4 Fie P(x0,y0)punct critic pentru funcţia f (x, y) şi 1, 2 minorii diagonali ai matricii Hessiene în acesta. Atunci: 1. Dacă 1 0 , 2 0 atunci P este punct de minim local pentru f (x, y); 2. Dacă 1 0 , 2 0 atunci P este punct de maxim local pentru f (x, y); 3. Dacă 1 0 , 2 0 sau 1 0 , 2 0 atunci P nu este punct de extrem local pentru f (x, y); 4. Dacă 1 0 sau 2 0 atunci nu putem preciza natura punctului critic P. Exemple : 10. Să se determine extremele funcţiei: f(x,y) x3 y2 3xy 15 131 Soluţie. Determinăm punctele de extrem local printre punctele critice adică printre soluţiile sistemului: f x 0, f y 0. Am văzut anterior că sunt două puncte critice (0,0), 9 27 , 4 8 . Fiecare poate să fie sau nu punct de extrem. 2f x2 6x, 2f x y 3, 2f y2 2 Matricea hessiană: H(x,y) 6x 3 3 2 O calculăm în fiecare punct critic: H(0,0) 0 3 3 2 şi 1 0, 2 9 Punctul (0,0) nu este punct de extrem. 9 27 H , 4 8 27 2 3 3 şi 1 2 27 , 2 2 18 Punctul 9 27 4 8 , este un punct de minim. 20. Să se determine extremele funcţiei: f(x1,x2,x3) x2 x2 x2 x1x3 1 2 3 Punctele staţionare (critice) sunt soluţii ale sistemului f x1 2x1 x3 0, f x2 2x2 0, f x3 2x3 x1 0 Obţinem un singur punct critic (0,0,0). Acesta este şi singurul posibil punct de extrem. 2f 2 1 x 2, 2f 2 2 x 2, 2f 2 3 x 2, 2f x1 y2 2, 2f x2 x1 0, 2f 2f x1 x3 x3 x1 1, 2f 2f x2 x3 x3 x2 0 132 Matricea hessiana devine: 2 0 1 H(x1,x2,x3) 0 2 0 1 0 2 şi are lanţul minorilor principali: 1 2, 2 4, 3 6 Rezultă că (0,0,0) este un punct de minim local. 133 Bibliografie 1. Allen, R.G.D., „Analiză matematică pentru economişti”, Ed. Şt., Bucureşti, 1971; 2. Dantzig, G., „Application et prologements de la programmation lineaire”, Dunod, Paris, 1966; 3. Diaconiţa, V., Manolachi, A., Rusu, G., „Matematici aplicate în economie”, Ed. Chemarea, Iaşi, 1993; 4. Diaconiţa, V., Rusu, G., „Optimizări liniare”, Ed. Sedcom Libris, Iaşi, 2001; 5. Diaconiţa, V., „Matematici aplicate în economie”, Ed. Paralela 45, Piteşti, 2003; 6. Diaconiţa, V., Rusu, G., Spînu, T.M., „Matematici aplicate în economie”, Ed. Sedcom Libris, Iaşi, 2004; 7. Drăgan, I., „Tehnici de bază în programarea liniară” , Ed. Tehnică, Bucureşti, 1976; 8. Mihoc, G., Ştefănescu, A., „Programarea matematică”, Ed. Did. şi Ped., Bucureşti, 1973; 9. Tamaş, V., „Matematică pentru studenţii economişti”, Ed. Junimea, Iaşi, 2001. 134 